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SUMMARY 

The.  thermodynamic  and  transport  piOpertieo  of  high.-temperd.ture  air 
are  found  in  closed  form  starting  from  approximate  partition  functions 
for  the  major  components  in  air  and  neglecting  all  minor  components. 

The  compressibility,  energy,  entropy,  the  specific  heats,  the  speed  of 
sound,  the  coefficients  of  viscosity  and  of  thermal  conductivity,  and 
the  Prandtl  numbers  for  air  are  tabulated  from  500°  to  15,000°  K  over  a 
range  of  pressure  from  0.0001  to  100  atmospheres.  The  enthalpy  cf  air 
and  the  mol  fractions  of  the  major  components  of  air  can  easily  be  found 
from  the  tabulated  values  for  compressibility  and  energy.  It  is  predicted 
that  the  i-randt!  HtunLei  for  Fully  io^iaed  air  will  become  small  compared 
to  unity,  the  order  of  0.01,  and  this  implies  that  boundary  layers  in 
such  flow  will  be  very  transparent  to  heat  flux. 


INTRODUCTION 


It  is  axiomatic  that  the  science  of  aerodynamics  must  be  based  oxi 
a  good  understanding  of  the  atmospheric  medium  through  which  vehicles 
are  to  fly.  Under  subsonic  flight  conditions,  air  may  be  treated  as  an 
gas  CampwCwu.  ai  rigid,  rotating  Hcitwijiic  molecules.  The  thermod^  — 
namic  properties  of  such  a  gas  are  well  known  and  they  are  accounted  for 
in  the  gas  flow  equations  by  the  familiar  ratio  <->f  specific  heats,  which 
in  this  case  is  a  constant.  Under  supersonic  flight  conditions,  air  may 
be  raised  to  temperatures  where  the  molecules  can  no  longer  be  treated 
as  simple,  rigid  rotators.  At  relatively  low  supersonic  speeds,  vibra¬ 
tional  energy  is  excited  arid  then  the  specific  heats  become  functions  of 
temperature.  However,  both  the  the rmo dynamic  and  transport  properties 
of  air  in  vibrational  excitation  can  be  predicted  with  fair  accuracy  by 
the  methods  01  quantum  statistics  and  Kinetic  tneory  (rei .  1) ,  ana  tne 
air- flow  relations  can  be  modified  accordingly.  Eggers  (ref.  2)  has 
calculated  the  effects  of  vibrational  energy  on  the  one -dimensional, 
inviscid  flow  of  diatomic  gases,  for  example. 

Further  changes  in  air  properties  may  occur  at  still  higher  flight 
velocity.  Flight  velocities  of  practical  interest  have  now  increased 
from  low  supersonic  speeds  to  the  escape  velocity,  37,000  feet  per  second. 
Vehicles  which  travel  at  these  hypervelocities  excite  the  air  to  such 
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high  temperatures  that  the  molecules  not  only  vibrate  but  may  dissociate 
Into  ate.,  i  and  even  ionize.  Under  these  conditions,  the  behavior  op  air 
deviates  widely  from  that  of  an  ideal  gas  and  the  thermodynamic  and  trans¬ 
port  properties  all  become  functions  of  pressure  as  well  as  of  tempera¬ 
ture.  It  is,  of  course,  essential  to  evaluate  these  functions  in  order 
to  calculate  the  pattern  of  air  flow  about  high-speed  vehicles,  the  vis¬ 
cous  and  pressure  forces  which  result,  and  the  heat  flux  which  occurs 
between  the  air  and  the  vehicle. 

The  equilibrium  thermodynamic  properties  of  a  gas  can  be  calculated 
with  great  confidence,  provided  the  energy  levels  of  the  gas  particles 
and  the  degeneracy  of  these  levels  are  known.  For  monatomic  and  diatomic 
gases  this  information  can  generally  be  deduced  from  spectroscopic  data 
with  such  accuracy  that  the  calculated  thermodynamic  properties  can  be 
trusted  to  very  high  temperatures,  even  where  experimental  confirmation 
is  lacking.  In  the  case  of  air,  however,  one  of  the  important  energy 
terms  was  not  known  with  confi dor -e  until  recently,  namely,  the  dissocia¬ 
tion  energy  of  molecular  nitrogt*;.  This  uncertainty  arose  because  the 
available  spectroscopic  data  were  consistent  with  two  different  models 
for  njtrogen  dissociation,  one  leading  to  a  dissociation  energy  of  7.37 
electron  volts  per  molecule  and  the  other  to  9*76  electron  volts  per 
molecule.  At  first,  the  lower  value  was  widely  accepted  as  the  most 
probable  one  (Herzberg,  ref.  3).  Krieger  and  White  (ref.  4)  and 
Hi r&chf elder  and  Jurli&s  (ref.  p)  ha/e  published  tablet,  of  thermodynamic 
properties  of  high- temperature  air  based  on  this  value.  Gaydon  (ref.  6) 
was  perhaps  one  of  the  first  advocates  of  the  view  that  the  higher  value 
was  the  correct  one.  Subsequently  a  number  of  experiments  were  performed 
which  confirmed  Gaydon’ s  opinion,  among  them  the  measurements  of  strong 
shock  waves  in  nitrogen  made  by  Christian,  Duff,  and  Yarger  (ref.  7)  and 
the  detonation  studies  made  by  Kistiakowsky,  Knight,  and  Malin  (ref.  8) . 
This  rendered  the  work  of  references  4  and  5  obsolete,  but  shortly  there¬ 
after  Gilmore  (ref.  9)  computed  the  chemical  composition,  energy,  entropy, 
compressibility,  and  pressure  of  air  as  functions  of  temperature  ana  den¬ 
sity  based  on  the  higher  value  for  the  dissociation  of  nitrogen.  Later, 
Hilsenrath , and  Beckett  (ref.  10)  published  a  similar  table  of  these  prop¬ 
erties,  but  in  much  smaller  increments  ox  temperature  and  density.  The 
calculations  in  both  of  these  references  (9  and  10)  are  highly  refined 
in  the  sense  that  they  not  only  account  for  the  major  components  of  air 
and  their  most  significant  energy  states,  but  they  also  take  into  account 
a  large  number  of  the  higher  energy  states  which  are  infrequently  excited, 
even  at  high  temperatures,  and  most  of  the  very  minor  chemical  components 
of  air  are  included.  Therefore  these  works  are  among  the  most  detailed 
estimates  for  the  therrAodynamic  properties  of  air  which  have  been  made. 

it  ia  likely  that,  these  will  bourne  since  the 

values  of  all  the  important  energy  levels  used  in  these  calculations  are 
now  quite  secure.  However,  it  is  desirable  to  have  approximate  expres¬ 
sions  for  these  properties  in  closed  form  which  can  he  solved  without 
iteration.  Such  solutions  would  be  particularly  valuable,  for  example, 
in  preparing  tables  to  be  used  with  the  method  of  characteristics  for 
calculating  the  flow  of  real  air. 
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In  contrast,  to  the  fairly  satisfactory  state  of  development  in 
regard  to  thensodynamlc  properties,  knowledge  of  the  transport  properties 
of  sir  at  high  temperatures  is  in  an  elementary  stage.  It  is  generally 
agreed  that  an  accurate  calculation  of  the  transport  properties  should 
he  based  on  the  rather  rigorous  theory  of  Chapman  and  Easkog  for  monatomic 
gases  (ref,  11).  The  extension  of  this  theory  given  by  Wang-Chang  and 
Uhlenbeck  (ref.  12)  would  be  used  for  molecular  gases  with  internal 
energy.  For  such  calculations  it  is  necessary  to  know  the  interaction 
potentials  which  exist  between  the  gas  particles  so  that  the  so-called 
collision  integrals  can  be  calculated.  H'irschfelder,  Curtiss,  and  Bird 
(ref,  13)  have  developed  the  theory  of  inte molecular  collisions  to  the 
stage  where  the  collision  integrals  and  the  transport  properties  of  inert, 
molecular  gases  can  be  calculated  with  good  accuracy.  However,  when  the 
air  dissociates,  as  it  does  at  high  temperatures,  the  atom-atom  and  atom- 
molecule  potentials  are  needed,  and  these  are  not  sufficiently  well  known 
to  calculate  the  collision  integrals.  At  the  present  time,  calculation 
of  the  atom-atom  and  the  atom-molecule  potentials  is  being  attempted  by 
quantum  mechanical  methods.  However,  it  may  be  som  time  before  these 
solutions  are  available,  and  even  when  they  are  compxeted  a  formidable 
obstacle  to  the  calculation  of  the  transport  properties  remains.  This 
occurs  because  two  atoms  may  approach  each  other  along  any  one  of  a  num¬ 
ber  of  potentials,  depending  upon  the  principal  quantum  number  of  the 
electrons  and  the  spatial  orientation  of  the  electron  spin  and  orbital 
momentum  vectors  at  the  time  of  collision.  These  potentials  are  quite 
different,  some  are  partly  attractive  and  others  are  entirely  repulsive. 

The  solutions  for  all  of  these-  potentials  are  implicit  in  the  quantum 
mechanical  description  of  the  problem.  It  is  hoped  that  these  potentials 
might  be  weighted  by  their  probabilities  to  yield  single  effective  poten- 
1 3 ale  - t  t-uch  Bpfeciffl,  for  If  the  at. me  which  interact  al.i.g  dif¬ 

ferent  potential  must  be  treated  as  a  separate  species  in  the  gas,  the 
calculations  of  the  transport  properties  will  be  laborious  indeed,  even 
with  electronic  computing.  Problems  of  this  type  are  further  compounded 
when  temperatures  are  considered  where  the  air  begins  to  ionize.  In  any 
event,  there  is  an  urgent  need  for  an  estimate  of  the  properties  of  air, 
due  to  the  demands  created  by  the  expanding  realm  of  very  high-speed 
flight.  Therefore,  an  engineering  approximation  for  the  transport  prop¬ 
erties  of 'Mgh-tiirp.i  raturt  air  wltjiM  be  valuable  in  the  interim  wh i la 
more  exact  solutions  are  being  prepared.  Even  after  more  precise  solu¬ 
tions  are  available,  an  approximation  giving  the  transport  properties  in 
closed  form  will  be  just  as  desirable  as  it  is  for  the  thermodynamic 
properties  of  air. 

* 

In  view  of  the  needs  outlined  above,  it  is  the  purpose  of  this  paper 
to  develop  approximate  expressions  for  the  properties  of  air  over  the 
range  of  temperatures  and  pressures  encountered  by  vehicles  traversaig 
the  atmosphere  at  speeds  up  to  escape  velocity.  The  principle  whicu  shall 
be  used  in  deriving  these  oxpressi^ns  is  that  thc.y  shall  be  made  as  -simple 
in  form  as  possible.  For  the  thermodynamic  properties  this  is  accomplished 
by  keeping  only  those  terns  which  are  necessary  to  yieV  a  final  result, 
within  a  few  percent  of  the  more  exact  solutions.  The  properties  of  air 
which  will  be  evaluated  are  the  compressibility  (i.e.,  the  correction  to 
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the  ideal  gas  equation  of  state),  the  energy,  enthalpy,  entropy,  specific 
heat  at  constant  pressure  and  at  constant  density,  the  speed  of  sound, 
the  viscosity,  thermal  conductivity,  and  the  Prandtl  number.  All  of 
these  properties  will  be  evaluated  for  equilibrium  conditions.  Because 
of  the  finite  reaction  rates,  these  values  will  doubtless  need  to  be  modi¬ 
fied  for  processes  which  involve  changes  In  sta^e  which  are  rapid  compared 
to  the  rate  of  approach  to  chemical  equilibrium.  Nonequilibri urn  effects 
will  probably  be  encountered  in  very  high-altitude  flight  because  the 
approach  to  equilibrium  is  slow  at  the  low  pressures  experienced  xhere. 
However,  there  is  experimental  evidence  that  equilibrium  is  essentially 
realized  in  flow  of  dissociating  air  under  conditions  encountered  in 
flight  at  moderate  altitudes  and  speeds  (ref.  14) .  In  addition, 

Hirschf elder  (ref.  15)  has  argued  that  heat  transfer  in  pure  conduction 
processes  will  correspond  to  equilibrium  values  if  the  reaction  rate  in 
one  direction  is  rapid,  Therefore,  the  thermodynamic  and  transport  prop¬ 
erties  of  air  which  are  based  on  equilibrium  conditions  should  apply 
directly  to  some  practical  problems  as  well  as  being  a  convenient  reference 
for  the  nonequilibrium  values. 

SYMBOLS 

speed  of  sound  (zero  frequency) 

stoichiometric  coefficients  for  components  and  Bp 
components  of  a  chemical  reaction 
Sutherland’s  constant  (eq.  (58)) 

specific  heat  per  mol  at  constant  density  for  component  i 
specific  heat  per  mol  at  constant,  pressure 

partial  specific  heat  per  mol  at  constant  pressure, 

£*i(Ci+l) 
i 

specific  heat  per  mol  at  constant  density 

dissociation  energy  per  molecule,  also  diffusion  coefficient 

_  %■  ~~ 

binary  diffusion  coefficient  for  molecules  of  type  i  into 

molecules  of  type  j 

base  of  natural  logarithms,  also  electron  charge 
electron 


a 

a-j  ,bt 

Ais®i 

C 

Ci 

Cp 


Cv 

D 

Dij 

e 

e~ 
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E 

Ei 

Eo 

Si 

Sn 


h 


H 

Hi 

I 

Io 

J 

k 

kn>^' 

kr 

KP 

Le! 


2n 

T1 

Mi 

M 

Mo 


energy  per  mol.  also  electric  field  strength 

energy  per  mol  of  component  i 

energy  per  mol  at  zero  absolute  temperature 

degeneracy  of  the  ith  state 

degeneracy  of  the  nth  electronic  state 

Planck*  s  constant 

enthalpy  per  mol 

enthalpy  per  mol  of  component  i 

molecular  moment  of  inertia,  also  ionization  energy  per 
molecule 

resonance  potential  for  ionization 
rotational  quantum  number 

Boltzmann  constant,  also  thermal  conductivity 

reference  coefficient  of  thermal  conductivity  (eq„  (77)) 

coefficient  of  thermal  conductivity  aue  to  molecular 
collisions 

coefficient  of  thermal  conductivity  due  to  chemical  reaction 
chemical  equilibrium  constant  for  concentration  units 
chemical  equilibrium  constant  for  pressure  units 

DpO* 

partial  Lewis  number ,  - £_ 

Mk* 

logarithm  to  the  base  e 

mass  of  a  gas  particle  =-  ® 

molecilar  weight  per  mol  of  component  i 

mean  molecular  weight  per  mol  of  a  gas  mixture 

molecular  weight  per  mol  for  undissociated  molecrJ.es 
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n 

n(Ai), 

n(Bi), 

N 

No 

N4 

N2 

NO 

0 

O4 

02 

P 

Po 

p(Ai), 

p(Bi), 

Pr 

Pr* 


Q 

Qt 

Qr 

Qv 

Qe 

Qc 


i 

w 

§ 

i 


vibrational.  quantum  number,  also  electronic  quantum  number, 
also  concentration  in  ao-1eu  pe"  unit  volume 

concent. *at ion  of  components  Ap,  Bq,  _ in  moles  per  unit 

volume 


nitroge'  atom,  also  atoms  in  general 

Avagadro  number,  vlecules  per  mol 

nitrogen  plus  ion,  niso  plus  ions  in  general 

nitrogen  m<.  lecule 

nitric  oxi 

oxygen  atom 

oxygen  plus  i  on 

oxygen  molecule 

pressure 

reference  pressure,  1  atmosphere 

partial  pressure  of  components  Aj_,  Bq,  ... 

iVq 

Prandtl  numBi  r,  -.i- 

!  I 


l 


f 

s 


partial  Prandtl  number. 


Mk1 


total,  partition  function 
translational  partition  function 
rotational  partition  function 
vibrational  partition  function 
electronic  partition  function 

total  partition  function  for  a  standard  state  of  unit 

concentration,  JL  Q 
’  RT 


me  A  TN  kiyO 
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Qp 


re 


R 

S 


Si 


So 

Sij  or  1 
S(i-j)  J 

T 

ui 

uo 

U 


xi 

x(Aj^)  . . . 
Z 


a 

3 


total  partition  function  for  a  standard  state  of  unit 
pressure,  pQ 

total  partition  functions  for  components  A^,  B±,  ... 
distance  between  atoms 

distance  between  atoms  where  the  potential  is  a  minimum 
universal  gas  constant,  energy  per  mol  deg 
entropy  per  mol 

entropy  per  mol  of  component  i  at  1  atmosphere  pressure 
collision  cross  section  for  undissociated  air  molecules 

collision  cross  section  for  particle  i  with  particle  j 
absolute  temperature 

mean  molecular  velocity  for  molecule  type  i 
mean  molecular  velocity  for  undissociated  air  molecules 
potential  energy  between  gas  particles 
mol  fraction 

mol  fraction  of  component  i 
mol  fraction  of  component  Aj_  ... 

pMn  Mq 

compressibility,  -  or  ~ 

pRT  M 

molecular  symmetry  number  (equal  2  for  homonuclear  diatomic 
molecules),  also  polarizability 

Morse  function  constant  (eq..  (6l)) 


7 

e 


Co 

ratio  of  specific  heats,  —a 

ef¬ 
fraction  of  molecules  which  are  dissociated  or  of  atoms 
which  sire  ionized 
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£i 

energy 

of 

the 

ith 

state 

£n 

energy 

of 

the 

nth 

electronic  state 

C 


dimensionless  distance  parar?  'VCf,  — 


1 


n 


7n 

Ai 

v 

Pi 

Po 

cr 


coefficient  of  viscosity- 

reference  coefficient  of  viscosity  (eq.  (67)) 
mean  free  path  for  molecule  type  i 
reference  mean  free  path  (eq.  (71)) 
vibrational  frequency 
density  of  molecule  type  i 
reference  density  (eq.  (70)) 
collision  diameter 


Subscripts 

p  partial  derivative  at  constant  pressure 

p  partial  derivative  at  constant  density¬ 
's  partial  derivative  at  constant  entropy 

ij j  indices  referring  to  molecules  type  i  and  j 

indices  referring  to  the  contribution  of  orans rational, 
rotational,  vibrational,  and  electronic  energy  modes, 
respectively 

<*• 

THERMODYNAMIC  PROPERTIES 


As  a  preliminary,  a  brief  review  of  some  of  the  results  of 
statistical  mechanics  will  be  given.  This  will  include  the  definitions 
of  the  partition  functions  and  will  summarise  those  relations  between 
these  functions  and  the  thermodynamic  properties  of  gases  which  will  be 
used  in  the  approximations  to  follow. 
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Partition  Functions 


All.  of  the  thermodynamic  properties  of  a  gas  may  he  calculated  from 
its  partition  function.  Consequently,  the  first  step  in  determining  the 
properties  of  air  is  to  calculate  the  partition  functions  for  the 
components  in  air.  The  partition  function  may  be  defined  as 


kT 


(1) 


where  is  the  energy  of  the  ith  state  of  the  gas  particle  end  g., 

is  the  degeneracy,  that  is,  the  number  of  states  of  the  particle  which 
have  this  same  energy  level.  The  energy  may  be  due  to  the  translational 
rotational,  or  vibrational  motion  of  the  particle,  or  to  the  motion  of 
the  electrons  within  the  particle.  The  temperatures  being  considered  in 
this  paper  are  in  all  cases  low  enough  that  the  excited  nuclear  energy 
states  may  be  disregarded.  The  usual  assumption  is  made  that  no  coupling 
exists  between  the  different  modes  of  energy.  Then  tne  partition  function 
may  be  expressed  as  the  product 


Q  =  Q+Qr-QyQe  (2a) 

The  factors  on  the  right  side  of  equation  (2a)  are,  respectively, 
the  partition  functions  associated  with  the  translational,  rotational, 
vibrational,  and  electronic  energy  levels  of  the  gas  particle.  Each  ’ 
factor  is  determined  independently  by  an  equation  of  the  same  form  as 
equation  (l).  By  the  methods  of  statistical  mechanics  it  is  found  that 
for  diatomic  molecules  these  factors  are : 


'  V  h*  J 


-nrjy 

p" 


CO 

Qr  = 

\ 

l 

(2J+l)e 

J=o 

00 

nhv 

<*v  = 

I 

e"  kT  = 

n=i 

00 

€n 

- 

I 

~  ~  kT 

8ne 

n=i 

b2J(  J+i) 

8Jt2IkT  _  Sj^IkT 


(2b) 


(2c) 


(2d) 


(2e) 
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The  nutations  used  above  are  common  ones  and  the  derivation  of  the 
equations  may  be  found  in  any  standard  Aext  on  statistical  mechanics, 
such  as  re fen 'nee  1 6  or  17.  For  monatomic  particles,  which  have  no  modes 
of  rotational  or  vibrational  energy,  the  rotational  and  vibrational  par¬ 
tition  functions  take  the  value  unity.  The  translational  and  electronic 
partition  functions  for  such  particles  take  the  same  form  as  equations  (?b) 
and  (2e),  respectively. 

Consistent  with  the  stated  purpose  of  this  paper,  only  those 
exponential  terms  are  included  in  the  electronic  partition  function 
(eq.  (2e) )  for  which  the  energy  levels,  en,  are  less  than  six  times  kT 
at  the  maximum  temperature  considered  (15,000°  K) .  Actually,  the  levels 
are  so  widely  split  in  this  range  that  the  closest  to  this  cutoff  is  the 
sixth  state  of  the  atomic  nitrogen  ion,  just  a  little  more  than  four 
times  kT  at  15,000°  K. 

Table  I  presents  the  atomic  and  molecular  constants  which  were  used 
in  calculating  the  partition  functions.  The  molecular  constants  for 
rotation,  vibration,  dissociation  energy,  and  electronic  energy  levels 
were  taken  from  Herzberg  (ref.  3).  The  9*76  electron-volt  value  for 
nitrogen  dissociation  is  used,  and  it  is  assumed  toat  the  rotational  and 
vibrational  constants  for  all  excited  electronic  states  are  the  same  as 
for  the  ground  state.  The  atomic  energy  levels  are  taken  from  Moore 
(ref.  18).  The  constants  have  been  rounded  off  to  the  nearest  0.5  per¬ 
cent,  and  the  second  and  third  electronic  energy  levels  of  atomic  nitro¬ 
gen  have  been  combined  since  they  lie  within  0.’5  percent  of  each  other. 

The  same  treatment  applies  to  the  second  and  third  electronic  levels  of 
the  atomic  oxygen  positive  ton. 

The  functions  which  are  to  be  used  directly  in  the  calculations  to 
follow  are  the  logarithms  of  the  partition  functions.  From  the  partition 
function  constants  (table  I)  and  equations  (2a)  through  (2e),  these 
functions  become 


to  Q(N2)  =  1  to  T  -  0.42  -  to(  1 


to  Q(02)  =  1  to  T  +  0.11  -  ln(  1 


11390 


to  (  3  +  2e  T  +  2e 


* 
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,  ZZB  _  326 

Zn  Q(0)  =  2  Zn  T  +  0.50  +  in  [  5  +  3e  ”T~  + 


22600  _  46600 

5e  ^  +  e  ^ 


-  2n  p 


(3c) 


27700 

in  Q(N)  -  “  in  T  +  O.30  +  Zn  I  4  +  XOe  T 


41500 

+  6e 


-  in  p 


•>+\  5 


38600  58200 


in  Q(0+)  Zn  T  t  O.50  +  Zn  [  4  +  10 


T  /- 
+  oe 


T 


(3d) 


-tap 


(3e) 


70.6  188.9 


Zn  Q(N+)  =  2  zn  T  +  0.30  +  in  (  X  +  3e  T  +  5e  T  + 

2  1 


22000  47000  67900 


5e 


T 


+  5e 


T 


in  p 


(3f) 


Zn  Q(e“)  =  ~  Zn  T  -  14„24  -  Zn  p  (3g) 


where  T  is  the  temperature  in  degrees  Kelvin  and  p  is  given  in 
atmospheres. 

It  will  he  surmised  that  those  components  of  air  which  are  not 
represented  hy  partition  functions  above  are  to  be  neglected  in  the 
approximations  which  are  to  follow.  The  absence  of  a  partition  function 
for  nitric  oxide  may  be  found  surprising,  but  it  will  be  seen  later  that, 
b  the  order  of  accuracy  being  considered  here,  the  formation  of  nitric 
oxide  may  be  neglected  in  computing  thermodynamic  properties  over  a  wide 
range  of  pressure  and  temperature  conditions,  including  those  conditions 
which  will  generally  be  encountered  in  high-speed  flight  through  the 
atmosphere.  This  occurs  because  NO  has  about  the  same  thermodynamic 
properties  as  an  average  for  nitrogen  and  oxygen  molecules,  and  the  NO 
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formation  does  not  great ly  Influence  the  equilibrium  mol  fractions  of 
molecules  and  of  atoms.  The  possible  influence  of  nitric  oxide  on  the 
transport  properties  of  air  will  be  considered  later. 


Energy,  Enthalpy,  and  Specific  Heats  for  the 
Components  of  Air 


According  to  statistical  mechanics,  the  energy  and  enthalpy  per  mol 
of  pure  gas  are  given  by  the  following  relations 


=  T 


=  T 


d  in  Qc 

(b) 

dT 

d  to  Qp 
dT 

(5) 

! 


The  quantities  Qc  and  Qy  are  the  partition  functions  for  the  standard 
states  of  unit  concentration  and  of  unit  pressure,  respectively.  These 
are  related  to  the  total  partition  function  by 


(6a) 

Qp  -  pQ 

(6b) 

and  they  are  functions  only  of  temperature  so  that  it  is  their  total 
derivatives  which  are  related  to  the  energy  and  enthalpy  as  given  in 
equations  (4)  and  (5)*  The  quantity  E0  is  a  constant  representing  the 
energy  of  the  gas  at  zero  absolute  temperature.  The  choice  of  this  level 
is  arbitrary,  but  by  convention  Eg  is  taken  as  zero  for  the  molecules 
of  nitrogen  and  oxygen.  Then  Eg  for  the  neutral  atoms  is  just  one  half 
the  energy  of  dissociation  per  mol  of  diatomic  molecules,  and  Eg  for 
the  xonized  atoms  is  the  sum  of  this  dissociation  energy  and  the  energy 
of  ionization.  With  all  the  ionization  energy  attributed  to  the  ionized 
atoms,  E0  for  the  electrons  must  be  taken  as  zero,  of  course. 

By  definition  the  specif,  heats  per  mol  of  pure  gas  are 


Cv  - 


fas') 

w, 


P 


dH 

Bt 


(7) 

(8) 


. 


■  I  x 
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*  ■■  '^«r 


Frum  equations  (2a),  (2b),  (2e)  (4), 
energy  and  enthalpy  per  mol  of  gas  due  to 
excitation  are  given  by 


and  (5)  it  is  seen  that  the 
translation  and  electronic 


(9) 


(10) 


and  the  specific  heats  given  by  equations  (j)  and  (8)  can  be  expressed 


(12) 


enthaWati0?SiV) '  (10)>  (U)’  6ive  that  part  uf  the  energy 

enthalpy  and  the  specific  heats  which  is  due  to  the  translational  Sd 

electronic  energy  for  either  atoms  or  molecules.  The  contributions  of 

rotations  and  vita***  ^  be  inalded“ 4  Se  “Lo^lar 

’  of  course.  According  to  equations  (2c)  and  (2d)  the  expression 


i-1 


w. 


=  1  + 


r+v 


(13) 


should  be  added  to  equations  (9)  and  (10) 
energy  and  enthalpy  for  diatomic  molecules 


in  order  to  obtain  the  total 
,  while 
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(iM 


should  he  added  to  equations  (ll)  and  (12)  to  get  the  total  specific 
heats  for  these  particles. 

The  entropy  of  a  pure  gas  is  related  to  its  partition  function  by 


S 

R 


(15) 


and  from  equations  (5)  Mid  (6b),  the  entropy  of  the  gas  at  unit  pressure 
may  be  expressed 


(16) 


The  calculated  values  of  the  partition  fi.  .tions  for  a  standard 
state  of  1,’ atmosphere  pressure,  the  energy,  and  the  specific  heat  at 
constant  density  for  the  major  components  of  air  are  given  in  table  II 
in  500°  K  increments  of  temperature.  The  concentration  standardized 
partition  functions,  the  enthalpy,  entropy,  and  the  specific  heat  at 
constant  pressure  are  not  listed  since  they  may  be  easily  found  from  the 
properties  tabulated  and  equations  (6),  (10),  (12),  and  (l6).  The  tabu¬ 
lated  properties  for  each  component  are  given  for  the  range  in  tempera¬ 
ture  in  which  the  component  exists  in  air,  at  pressures  between  0.0001 
and  100  atmospheres,  up  to  the  maximum  temperature  considered,  15,000  ;  K. 
We  shall  now  examine  how  these  values  may  be  used  in  computing  the 
thermodynamic  properties  of  air. 


Equilibrium  Constants 


In  order  to  determine  the  equilibrium  mol  fractions  for  the  components 
of  air,  it  will  be  necessary  to  calculate  the  equilibrium  constants  for 
the  chemical  reactions  which  occur.  These  chemical  reactions  may  be 
expressed  in  the  general  form 


(IT) 


where  the  A±  are  the  reactants,  the  Bi  the  products,  and  &±  and  bq 
are  their  respective  stoichiometric  coefficients.  The  pressure  equilib¬ 
rium  constant  for  this  reaction  is  defined  in  terms  of  the  partial 
pressures 
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•  „nPbl(Bj) 
p  nPRi(Ai) 


(18) 


and  it  is  related  to  the  partition  functions  by  (see  ref.  16) 

to  Kp  -  -  +  £bi  hi  QpCBiJ-^ai  Zn  Qp(Ai)  (l 9) 

where 

BofBi)-  ^^^(Ap) 


is  the  zero  point  energy  of  the  products  less  the  zero  point  energy  of 
the  reactants,  both  referred  to  their  standard  states.  The  reactions 
considered  here  are  the  dissociation  of  molecular  oxygen  and  of  molecular 
nitrogen  and  the  ionization  of  atomic  oxygen  and  of  atomic  nitrogen . 
According  to  equation  (19)  and  the  partition  function  constants  given  in 
table  I,  the  equilibrium  constants  for  these  reactions  are 


W  -*2  0) 

Zn  Kp(N2  2N) 

In  Kt)(0  -■>  0++e“) 

c 

Zn  Kp(N  I^+e") 


—  +  ^  hi  *  to  'tp('v-s)  (20a) 


-3^°°  +  2  Zn  Qp(N)  -  Zn  Qp(N2)  (20b) 


^8?0Q°  +  Zr 1  4p(0+)  +  Zn  Qp(e-)  -  Zn  Qp(0) 

T  (20c) 

-%8°°  +  to  Qp(w!)  +  Zn  Qp(e-)  -  Zn  Qp(N) 

T  (20d) 


The  concentration  equilibrium  constant  is  defined  by 


nnbi(Bj) 

nnai(Ai) 


(21) 


where  n(Ap)  and  n(Bp)  are,  respectively,  the  concentrations  of  the 
chemical  reactants  and  products .  This  quantity  will  also  be  needed  for 
subsequent  calculations ,  and  it  is  obtained  by  replacing  the  pressure 
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standardized,  partition  functions,  Qp,  with  the  corresponding  concentration 
standardized  partition  functions,  Op  (eq.  (19)).  From  equation  (b)  it 
is  seen  that 


\ 


Kc  =  Kp(RT) 


£ai-£hi 


(22) 


The  logar i third c  derivatives  of  the  equilibrium  constants  will  also 
be  required  later.  From  equations  (5),  (10),  and  (19)  these  become 


T 


d  Zn  Kc 


(23) 


T  d  KP  -  T  d  Kc 
dT  dT 


(24) 


The  equilibrium  constants  and  their  logarithmic  derivatives  for  the 
reactions  represented  by  equations  (20a)  through  (20d)  are  listed  as 
functions  of  temperature  in  table  III.  These  quantities  will  now  be 
used  in  calculating  the  component  mol  fractions  and  their  derivatives. 


Calculation  of  the  Equilibrium  Mol  Fractions 
and  Their  Derivatives 


The  possibility  that  approximate  solutions  in  closed  form  could  be 
obtained  for  the  properties  of  air  suggests  itself  upon  examination  ol 
the  results  of  Gilmore  (ref.  9) •  His  tables  of  the  composition  of  air 
show  that  there  are  four  chemical  reactions  of  major  importance.  These 
are  the  dissociation  of  molecular  oxygen  and  of  molecular  nitrogen,  and 
the  ionization  of  atomic  oxygen  and  of  atomic  nitrogen. 


02  2  0 

(25a) 

N2  -as*  2N 

(25b) 

* 

0  — >  0+  +  e” 

(25c) 

N  N+  +  e~ 

(25d) 

With  one  exception,  all  other  reactions  which  occur  yield  component 
concentrations  which  are  the  order  of  0.1  percent,  or  less,  of  the  major 
components  given  by  the  reactions  above.  The  exception  is  the  formation 

NACA  TN  4150 


17 


of  nitric  oxide,  NO,  which  i*  sea  level  density  may  become  as  much  as 
10  percent  of  the  air  around  c'000°  K  (see  ref.  9).  However,  even  this 
much  nitric  oxide  does  not  strongly  influence  the  resulting  thermodynamic 
properties  of  air,  and  at  densit  ^  less  than  0.01  normal  sea  level 
density,  where  the  NO  is  less  t nan  1  percent,  the  effects  are  very 
small . 

Two  distinctive  features  of  the  chemical  reactions  given  above  are 
observable  from  Gilmore's  results.  The  first  is  that  at  all  pressures 
the  dissociation  of  oxygen  is  essentially  complete  before  the  dissocia¬ 
tion  of  nitrogen  begins.  This  means  that  these  two  reactions  can  be 
treated  independently  for  the  purposes  of  cur  approximation.  The  stcond 
feature  is  that  nitrogen  and  oxygen  atoms  ionize  at  about  the  same  tem¬ 
perature  and  with  about  the  same  energy  changes.  (Note,  e.g.,  the  approx¬ 
imate  ejuaUt.gr  nf  equilibrium  cvvs+snts  ter  tSie  two  ionisation 
reactions,  table  XII.)  Consequently,  it  is  possible  to  assume  that  once 
air  is  completely  dissociated,  all  atoms  constitute  a  single  species 
which  has  the  population  weighted  average  properties  of  the  nitrogen  and 
oxygen  atoms.  Then  the  problem  conveniently  divides  itself  into  three 
cases,  and  in  each  case  only  three  components  of  the  air  need  to  be  con¬ 
sidered  simultaneously.  These  three  cases  are  described  as  follows, 
where  the  ratio  of  nitrogen  to  oxygen  has  been  taken  as  4  to  1. 

Case  I.-  At  relatively  low  temperatures,  the  three  components  of 
air  which  may  exist  simultaneously  are  molecular  nitrogen,  molecular 
oxygen,  and  atomic  oxygen.  If  e  is  the  fraction  of  the  initial  number 
of  moles  per  unit  volume  which  have  become  dissociated,  then  the  partial 
pressures  for  these  three  components  may  be  expressed 


p(Na) 

=  x(N2)p  =  M  p 

1+e 

(26a) 

p(o2) 

=  x(0p)p  =  °;2-e  p 

1+e 

(26b) 

p(0) 

=  x(0)p  =  —  p 

1+e 

(260) 

where  the  x(AjJ  are  the  mol  fractions  of  component  Aj_.  It  is  assumed 
that  the  gas  components  each  behave  like  an  ideal  gas.  Then  the  equation 
of  state  is 


p  _  ZRT 
P  “  Mq 


(27) 


where  the  compressibility  Z  is  1  +  e  in  this  case.  This  quantity  Z 
is  the  total  number  of  moles  per  initial  mol  of  undissociated  air.  It 
js  also  equal  t:  the  rati  j  jf  the  initial,  molecular  of 

undissociated  air  to  the  mean  molecular  weight,  Mq/S. 


NAC/1  ”f  UVjO 


IS 


The  equilibrium  constant  for  the  oxygen  dissociation  reaction  is 


Kp(02  -*20)  « 

p(02) 


4cgp 

(l+e) (0.2-e) 


(28) 


Then  the  fraction  €  is  found  by  solving  the  quadratic  equation  (28) 


(29) 


The  condition  of  no  dissociation  which  occurs  at  low  temperatures  "is, 
of  course,  just  the  limit  of  case  I  where  e  is  zero. 


Case  II.-  At  intermediate  temperatures  the  oxygen  is  completely 
dissociated  and  the  nitrogen  dissociation  commences.  For  this  case,  e 
will  be  defined  as  the  fraction  of  the  initial  air  which  dissociates  into 
nitrogen  atoms.  Then  tne  compressibility  Z  is  1.2  +  e  for  this  case, 
and  the  partial  pressures  of  the  components  are 


p(n2) 

p(N) 

p(o) 


=  x(N2)p  =  —  P 
1 .£+e 

=  x(N)p  =  p 

1.2+e 

=  x(c)p  =  p 

1.2+e 


(30a) 

(30b) 

(30c) 


The  equilibrium  instant  for  the  nitrogen  dissociation  reaction  is 

k_(n2  2N)  =  _ 

p  p(N2)  (1.2+e) (0.8-e) 


(31) 


whence  e  is  given  by 
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At  high  temperatures,  the  dissociation  of  nitrogen  is 
fills’'  complete ,  and  the  ionization,  of  the  atoms  begins.  For  this  case 
’.e  assume  that  the  atoms  are  a  single  homogeneous  species  symbolized 
by  N.  The  term  e  will  now  be  the  fraction  of  the  atoms  which  have 
been  ionized,  and  then  +%»e  partial  pressures  become 


p(N)  =  x(N)p  =  —  p 

1+e 

(33a) 

p(N+)  =  x(N+)p  =  JL-  P 

1+e 

(33b) 

p(e“)  =  x(e")p  =  p 

1+e 

(33c) 

Note  that  the  compressibility  Z  in  this  case  is  2(l+e).  The  ionization 
reaction  equilibrium  constant  is 


Kp(N  -*•  N^+e-)  =  P(Ij+)p(e")  =  _egP 

p(N)  1-e2 

and  then  e  becomes 


The  derivatives  of  the  mol  fraction  compressibility  products,  Zx±, 
will  be  needed  in  addition  to  the  mol  fractions.  For  case  I,  from 
equations  (28),  (26a),  (26b),  and  (26c)  it  is  seen  that 


(34) 

|  t 
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where  the  subscript  p  refers  to  the  fact  that  the  partial  derivatives 
are  taken  at  constant  pressure.  In  order  to  find  these  same  partial 
derivatives  at  constant  density,  the  equilibrium  constant  is  expressed 
in  concentration  units 


whence 


Kc(02  20) 


4egn _ 4e2  p 

(l+e)(0.2-e)  8:5  0.2-e  Mq 


d  Zn  Kc  /2  1  ^  AA 

dT  “  \€  +  0.2-e J  \bT/0 


(to) 

m 


and  the  mol  fraction  derivatives  at  constant  density  take  the  same  form 
as  equations  (37),  (38),  and  (39)  except  that  the  quantity  (de/dT)p  is 

replaced  by  the  quantity  (de/dT)p  found  in  equation  (41).  For  case  II 
a  similar  set  of  equations  for  the  derivatives  occurs 


d  Kp  _  (2  _  1  1 

dT  ”  V e  1.2+e  +  0.8-e, 

d  to  Kc  _  (2  1  \  fbe] 

dT  V€  0.8-e/  W/p 

cMn)  =  _2  azx(w2)  =  2  de 

dT  dT  dT 

dM°)  =  0 
St 


(42a) 

(42b) 

(43) 

(44) 


Finally,  for  case  III,  the  equations  are 


1  ln  Kc  ..('1+  Jj\/V\ 

dT  \€  1-e/  \§T/p 


dZx^'fl*)  _  dZx(e“)  _  _  dZx(N)  _  dg 
St  St  dT  dT 


(45a) 

(451) 

(46) 


#/here,  as  before,  the  psnsxal  u.ez'xvahi'vcs  are  given  for  constant  pressure 
or  for  constant  density  depending  on  the  corresponding  partial  derivative 
of  e  which  is  used.  With  these  quantities  in  hand  we  are  in  position 
to  calculate  the  energy,  specific  heat,  entropy,  and  the  speed  of  sound 
for  air. 
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Energy,  Specific  Heat,  Entropy,  and  Speed  of  Sound 
for  Air  in  Equilibrium 


The  energy  per  mol  of  air  is  simply  the  sum 


-I 


xi% 


(4f) 


where  the  Ei  are  the  energy  per  mol  for  component  i.  Generally  in 
aerodynamic  problems  the  energy  per  ft**  mass  of  gas  is  needed  rather 

ST™?6  per,mol‘  The  ““«*  of  in  a  mol  of  undissociated  air 

'rh.L  k  i}1  be,us1ed  as  the  constant  reference  quantity  for  this  purpose. 
Ihen,  in  dimensionless  form,  the  energy  per  mol  of  initially  undissociated 

air  xS 


s.zy** 

RT  L>  RT 
i 


(48) 


while  the  dimensionless  enthalpy  per  initial  mol  of  air  becomes 


ZH  =  ZE 
RT  RT  + 


(49) 


The  compressibility,  Z,  and  the  dimensionless  energy,  ZE/RT  which 

*£?  TewCfC^ir°"  thS  equations ,  are  Sted  in’ 

tables  17(a)  and  IV(b>  and  are  graphed  as  functions  of  temperature  in 

figures  1  and  2.  The  enthalpy  i5  easily  obtained  from  these  values  via 

nf^pf011  ^  *o  At,b^is  point  it  is  desirable  to  compare  the  calculations 
f  references  9  ana  10  with  the  results  obtained  above.  It  is  not  possi- 
ble  to  compare  the  numerical  results  directly  since  both  references  9 

?!  properties  of  air  as  Actions  of  temperature  at  con- 
stant  density  rather  than  at  constant  pressure  as  is  done  here.  However, 
the  results  uf  these  two  references  are  the  same,  and  reference  10  tabu¬ 
lates  the  properties  in  close  enough  intervals  so  that  constant  pressure 
.  lues  can  be  obtained  fairly  accurately  by  interpolation.  The  values 
or  compressibility  and  energy  interpolated  from  reference  10  agree  with 
the  approximate  solutions  within  5  percent  in  all  cases,  and  generally 
within  2  percent,  so  that  on  the  scale  of  figures  1  and  2  the  differences 

the  tT  discernibbe‘  In  vlew  °f  this  agreement,  it  is  concluded  that 
the  solutions  presented  here  a„re  useful  Arm™  vjhsst ■»  n  +t — 

S?  Pr0rrtlr  f  alr-  sh^l  pt^eeTtolse’thosl^a^rox- 

lff10'18  to  calculate  the  entropy,  specific  heats,  and  speed  of  sound  In 

al  x  i 


The  entropy  per  initial  mol  of  air  is  simply  obtained  from  the 
entropies  of  the  components  of  air  through  the  relation 
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where  pn  is  the  reference  pressure  for  the  standard  state,  in  this  case, 
1  atmosphere.  The  entropy  values  are  listed  in  table  IV(c)  and  are  shown 
graphically  in  figure  3, 

It  is  necessary  to  take  the  derivatives  of  the  products  ZE  and  ZH 
for  specific  heats  if  the  latter  are  to  have  their  usual  meanings,  that 
is  the  change  in  energy  which  occurs  in  a  fixed  mass  of  gas  per  degree 
temperature  change.  Thus  the  constant  density  specific  heat  per  initial 
mol  of  air  is  given  by 


where  Ci  is  the  derivative  of  energy  for  component  i,  that  is,  dEi/dT. 
The  corresponding  specific  heat  for  constant  pressure  is  given  by 


The  specific  heats  calculated  from  equations  (51)  (52)  are 

listed  in  tables  IV(d)  and  IV(e).  It  may  be  noted  that  when  chemical 
reactions  occur,  the  difference  between  the  specific  heats  per  mol  is 
not  equal  to  the  gas  constant,  as  in  an  ideal  gas,  but  is  given  by 


Both  of  the  specific  heat  functions  have  the  same  general  features. 
These  features  are  illustrated  in  figure  4  where  the  specific  heat  at 
constant  density  is  given  as  a  function  of  temperature  for  pressures 
of  0.0001  0.01,  and  1  atmosphere.  At  low  temperature  Cv  increases 

from  5/2  R  to  7/2  R  as  the  vibrational  modes  of  energy  become  excited. 
Then,  with  increasing  temperature,  the  specific  heat  goes  through  three 
distinct  maxima  where  the  chemical  components  change  most  rapidly  with 
temperature;  the  first  maximum  is  due  to  the  oxygen  dissociation  reaction, 
the  second  to  the  nitrogen  dissociation,  and  the  third  to  the  ionisation 
reactions.  As  pressure  decreases,  these  maxima  increase  in  sharpness  and 
in  magnitude  and  they  shift  to  lower  temperatures.  It  has  been  pointed 
out  (ref.  16)  that  the  values  determined  from  the  partition  functions 
cannot  De  highly  accurate  because  the  second  derivatives  of  the  partition 
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functions  are  required,  and  the  error  in  neglecting  higher  order  terms 
is  magnified  in  the  process.  However,  the  curves  of  figure  4  would 
merely  be  shifted  slightly  in  temperature  and  have  somewhat  different 
peak  values  if  the  calculations  were  more  exact.  Certainly  the  qualita¬ 
tive  features  of  the  specific  heats  should  be  correctly  g.iv  ;x. 

The  specific  heat  values  determined  above  enable  us  to  calculate 
the  speed  of  sound  in  air.  This  speed  of  sound  will  be  defined  by 


(54) 


that  is,  the  limiting  value  as  the  frequency, of  sound  approaches  zero. 
The  partial  derivative  given  in  equation  (54)  will  not  be  calculated 
directly,  since  it  is  not  convenient  to  treat  the  entropy  as  an  inde¬ 
pendent  variable  which  can  be  held  constant  (see  eq.  (50)).  However, 
equation  (54)  may  be  transformed  with  perfect  generality  into  the  form 


where  7  is  the  ratio  of  the  specific  heats,  Cp/Cv.  Equation  (55)  is, 
in  turn,  equivalent  to 


which,  from  the  equation  of  state  (eq.  (27)),  may  be  expressed  in  terras 
_,f  vstrial  LcS  whlrvh  have  already  been  calculated 


The  dimensionless  speed  of  sound  parameter,  a2p/p,  is  listed  in  table  IV\f) 
and  is  plotted  as  a  function  of  temperature  in  figure  5.  The  second  term 
on  the  right  side  of  equation  (57)  is  generally  near  unity,  so  that  fig¬ 
ure  5  is  also  indicative  of  the  variation  in  7  with  temperature. 

The  range  of  pressure  and  temperature  in  which  the  dissociation  of 
oxygen,  the  dissociation  of  nitrogen,  and  the  ionisation  cf  afccmi  are  the 
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important  chemical  reactions  is  indicated  by  dividing  table  IV  into  the 
regions  case  I,  case  II,  and  case  III,  respectively.  Figure  6  presents 
the  same  information  in  a  graphical  form  from  which  it  is  convenient  tc 
judge  the  chemical  state  of  air  for  a  given  pressure  and  temperature. 

It  will  be  noted  that  boundaries  between  the  chemical  reactions  corre¬ 
spond  to  end  points  in  the  charp.cteristic  cycles  that  occur  in  the 
thermodynamic  parameters  graphed  in  figures  1  through  5* 


Aerodynamic  Considerations 


The  thermodynamic  properties  obtained  at  this  point  are  those 
required  to  perform  calculations  of  inviscid  air-flow  problems  by  the 
method  of  characteristics.  These  properties  are  given  for  a  range  of 
temperature  from  500°  to  15,000°  K  and  of  pressure  from  0.0001  to 
100  atmospheres.  It  is  of  interest  now  to  examine  the  altitude  and 
velocity  at  which  these  conditions  will  occur  in  flight.  A  grid  of  the 
pressure  and  temperature  at  the  stagnation  point  of  a  body  in  flight  is 
shown  in  figure  7  as  a  function  of  flight  altitude  and  velocity.  The 
stagnation  enthalpy  per  unit  mass  was  simply  taken  as  one  half  the  veloc¬ 
ity  squared,  and  the  stagnation  pressure  was  related  to  the  static  pres¬ 
sure  (and  thus  to  aloitude)  with  the  results  of  Feldman  (ref.  19),  who 
computed  the  pressure  ratio  developed  across  normal  shock  waves  in  air 
at  various  altitudes.  Generally  lower  temperatures  and  pressures  will 
be  attained  at  regions  other  than  the  stagnation  region,  so  the  range  of 
these  variables  will  be  adequate  for  such  cases  also.  At  very  low  pres¬ 
sures,  such  as  occur  at  high  altitudes  or  in  expanded  flow  regions,  aero¬ 
dynamic  effects  may  often  be  disregarded.  However,  if  necessary,  the 
approximations  of  this  report  may  be  extended  to  pressures  below  0.0001 
atmosphere.  It  is  not  advisable  to  extend  them  to  pressures  much  greater 
than  100  atmospheres,  since  there  occurs  an  increasing  overlap  of  the 
chemical  reactions  in  air  as  pressure  increases  and,  in  addition,  the 
nitric  oxide  becomes  such  a  large  component  of  the  air  that  it  must  be 
considered.  As  a  result  of  this  overlap  in  the  chemical  reactions,  the 
calculated  properties  of  air  show  irregularities  in  the  region  of  transi¬ 
tion  from  one  reaction  to  another.  Slight  irregularities  can  be  observed 
at  the  transition  regions  for  10  and  100  atmosphere  pressures  (figs.  1 
through  4),  for  example-  However,  the  deviations  are  s.’iall  at  these  pres¬ 
sures  and  it  can  be  seen  from  figure  7  that  the  thermodynamic  properties 
of  air  can  be  closely  approximated  in  closed  form  over  the  range  of  con¬ 
ditions  of  current  interest  in  aerodynamics.  Some  of  the  results  of  these 
approximations  will  next  be  used  in  estimating  the  transport  properties 
of  high-temTJerature  air. 
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TRANSPORT  PROPERTIES 
Collision  Cross  Sections 


Consider  first  some  qualitative  aspects  of  the  collisions  between  gas 
particles.  The  particle  trajectories  are  influenced  by  a  potential  U 
which  is  negligible  at  long  range,  which  may  be  attractive  or  repulsive 
at  intermediate  range,  but  which  always  becomes  strongly  repulsive  at 
very  short  range.  A  particle  with  kinetic  energy  kT  will  not  be  greatly 
deflected  if  it  passes  only  through  that  part  of  the  potential  field 
where  (u (  «  kT.  On  the  other  hand,  the  particle  will  he  deflected  appre¬ 
ciably  If  It  parses  vfrew  (  »  kT.  The  direction  of  the  deflection  is 
unimportant  so  far  as  transport  properties  are  concerned,  for  it  is  the 
absolute  value  of  the  deflection  which  determines  the  change  in  mass, 
momentum,  and  energy  fluxes  caused  by  the  collision.  To  a  first  approxi¬ 
mation,  the  absolute  value  of  this  deflection  is  independent  of  the  sign 
of  the  potential  and  the  effective  collision  diameter  o'  is  the  order 
of  the  largest  distance  where  U  =  ±kT.  The  effective  collision  cross 
section  S  will  be  defined  as  n a2. 

In  the  rigorous  treatment  of  the  transport  properties  of  gases,  the 
effective  collision  cross  section  S  is  found  to  be  an  integral  function 
of  the^def J^ct ion  angles  produced  by  a  collision,  and  this  integral,  is  a 
functTonof  the  relative  velocity  of  the  colliding  particles.  The 
so-called  "collision  integral"  is  a  function  of  temperature  only  which, 
is  S  times  a  velocity  function  integrated  over  al  1  velocities.  Thus 
the  collision  integral  may  be  thought  of  as  an  average  effective  colli¬ 
sion  cross  section,  and  the  transport  coefficients  can  be  related  directly 
to  these  integrals.  However,  not  all  the  interparticle  potentials  have 
been  developed  which  are  needed  to  calculate  the  collision  integrals  for 
air.  In  the  present  paper  then,  plausible  estimates  of  the  effective 
collision  cross  sections  will  be  used  to  determine  the  mean  fre6:j3aths 
for  hard  elastic  spheres,  and  for  such  particles  the  transport 
cients  can  he  related  to  these  mean  free  paths  (ref.  1).  The  effects  of 
the  interaction  potentials  will  be  taken  into  account  by  letting  the 
spherical  cross  sections  be  appropriate  functions  of  temperature. 

The  effective  cross  sections  for  collisions  between  diatomic 
molecules  can  be  obtained  quite  accurately  by  the  collision  integral 
method.  However,  at  high  temperatures,  the  very  steep  repulsive  portion 
of  the  intermole cular  potential  is  penetrated  so  that  the  molecules  do 
behave  essentially  like  hard  spheres.  Consistent  with  the  approximations 
which  are  to  follow,  it  will  be  sufficient  to  use  the  Sutherland  formula 
fcr  the  molecular  cross  section  So 

+  |  (58) 

^00 

where  C  is  112°  K  and  S^,,  the  molecular  cross  section  at  infinite 
temperature,  is  3.l4xl0~l5cm2. 
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For  atomic  collisions,  the  picture  is  complicated  by  the  fact  that 
two  adorns  may  approach  each  other  along  any  one  of  a  number  of  potentials. 

For  example,  the  potentials  between  two  nor¬ 
mal  nitrogen  atoms  are  shown  qualitatively 
in  sketch  (a)  (see  ref.  3)*  The  lowest  lying 
of  these  potentials,  designated  UQ,  has  the 
lowest  total  electron  spin  and  it  is  the  one 
normally  responsible  for  the  vibrational 
energy  levels  observed  in  the  stable  molec¬ 
ular  state.  Therefore  UQ  can  be  expressed 
quantitatively  from  experimental  spectro¬ 
scopic  data.  Unfortunately,  the  higher  lying 
potentials  for  the  atoms  in  air  are  not  known 
quantitatively  at  present,  so  we  are  reduced 
to  estimating  an  average  collision  cross 
section  for  all  of  the  potentials  by  means 
of  the  known  lowest  lying  potential  U0. 

For  this  purpose  :Lt  is  assumed  that  the 
collision  diameters  a  are  given  by 


Sketch  (a) 


U0(c)  =  -kT  (59) 

and  these  diameters  will  be  used  to  evaluate  the  coefficients  of  momentum 
and  energy  transfer.  A  somewhat  deeper  penetration  of  the  potential  is 
normally  required  for  a  collision  to  affect  the  particle  flux,  so  that 
the  diameters  a*  which  will  be  used  to  evaluate  the  diffusion 
coefficients  are  assumed  to  be  given  by 


U0(a’)=-2kT  (60) 

It  may  be  pointed  out  that  Hirschf elder  and  Eliason  (ref.  20)  have 
examined  the  relation  between  values  of  the  transport  coefficients  given 
by  the  hard  sphere  model  and  by  the  more  rigorous  collision  integral 
method.  They  find  that  U(cr)  and  U(cr')  are  about -0.6kT  for  a  wide  vari¬ 
ety  of  attractive  potentials,  that  U(cr)  is  about  0.9kT  and  U(u* )  is 
about  1.6kT  for  a  similar  variety  of  repulsive  potentials.  If  all  the 
interparticle  potentials  Mere  known,  it  would  be  simple  to  use  these 
criteria  to  obtain  a  weighted  average  collision  diameter.  All  the  poten¬ 
tials  are  known  for  two  normal  hydrogen  atoms  (ref.  21)  and  it  is  found 
for  this  case  that  equations  (59)  and  (60)  yield  values  for  a  and  cr ’ 
which  agree  with  the  weighted  average  collision  diameters  within  8  per¬ 
cent  over  the  range  of  temperatures  from  1000°  to  15,000°  K.  Of  course, 
there  is  no  assurance  that  these  same  relations  will  hold  as  closely  for 
collisions  between  the  atoms  in  air.  In  fact  the  average  collision  diam¬ 
eters  for  normal  oxygen  and  nitrogen  atoms  will  probably  be  overestimated 
by  equations  (59)  and  (60),  since  the  shallow  intermediate  potentials 
(such  as  Ui  and  U2,  sketch  (a))  must  be  considered  for  these  atoms, 
whereas  they  do  not  occur  for  hydrogen.  The  effect  of  these  intermediate 
potentials  will  be  partly  compensated  for  by  the  fact  that  some  of  the 
atoms  will  be  in  excited  electronic  states  which  have  collision  diameters 
the  order  of  three  times  larger  than  the  normal  atoms  (ref.  20).  The 
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fraction  atoms  in  excited  states  is  small  over  most  of  the  temperature 
and  pressure  range  considered  so  that  collisions  between  excited  parti¬ 
cles  are  relatively  rare.  The  encounters  between  an  excited  and  a  normal, 
atom  are  the  ones  which  significantly  influence  the  mean  free  paths,  and 
the  cross  sections  for  such  collisions  are  about  four  times  normal.  The 
fraction  of  each  constituent  in  air  which  is  excited  is  listed  as  a 
function  of  temperature  in  the  following  table; 


— 1 

°K 

Fraction  of 
oxygen  atoms 
excited 

Fraction  of 
nitrogen  atoms 
excited 

Fraction  of 
ionized  atoms 
excited 

4,000 

0.002 

0.002 

6,000 

-013 

.025 

0.012 

8,000 

.032 

.031 

•  032 

10,000 

.055 

•  153 

.060 

12,000 

14,000 

.079 

.228 

.089 

.120 

ihe  first  twu  excited  states  of  atomic  oxygen  and  of  the  atomic  nitrogen 
plus  ion  are  very  close  to  the  ground  state  (see  table  I)  and  so  they 
have  not  been  counted  as  excited  states.  It  can  be  seen  that  the  frac- 

\,f  excited  py.rtido.cs  in  air  will  gti.ciallj  be  le&s  than  lu  percent 
and  in  the  very  worst  case  considered  about  20  percent  of  the  atoms  will 
be  excited.  This  occurs  at  100  atmospheres  pressure  where  the  ionization 
reaction  is  repressed  until  temperatures  beyond  12,000°  K  are  reached  and 
23  percent  of  the  nitrogen  and  8  percent  of  the  oxygen  atoms  are  excited. 
For  a  pure  gas  in  which  excited  particles  have  three  times  normal  colli¬ 
sion  diameters,  the  average  mean  free  path  is  decreased  28  percent  when 
10  percent  of  the  particles  are  excited,  and  is  decreased  45  percent 
when  20  percent  of  the  particles  are  excited.  Thus  the  error  introduced 
by  neglecting  the  long-range  forces  of  the  excited  atoms  is  probably  the 
same  order  of  magnitude  and  opposite  in  effect  to  the  uncertainties  caused 
by  neglecting  the  shortened  range  of  the  shallow,  intermediate  collision 
potentials . 

The  lowest  lying  atom-atom  potentials  may  be  approximated  by  the 
Mo  rise  function 


where  D  is  the  dissociation  energy  and  (;  is  the  dimensionless  distance 
parameter  (r/re-l) ;  re  ’is  the  interatomic  distance  at  which  UQ  is  a 
minimum  (sketch  (a));  3  is  a  constant  related  to  the  vibrational  fre¬ 
quency  v  of  the  stable  molecule  and  is  given  by  v(  1^1/20) 1/2.  Although 
the  Morse  function  is  not  very  accurate  at  long  ranges,  it  will  be  needed 
only  at  rather (high  temperatures  where  the  potential  is  fairly  well 
described  for  |U0|  >  kT.  The  Morse  function  constants  were  taken  from 
Herzberg  (ref.  3)  and  these  are  listed  in  the  following  table: 


Atomic 

particles 

X’  * 

re, 

angstroms  j 

0-0 

59,000  3*24 

1.207 

N-N 

113,200  2.96 

1.094 

N-0 

75,400  3*18 

1.151 

0-0+ 

75,200  3.18 

1.123 

N-N+ 

101,200  2.94 

1.126 

From  equations  (59)  and  (6l),  the  momentum  and  energy  transfer 
collision  diameter  a  is  given  by 


(62a) 


and  by  equations  (60)  and  (6l)  the  diffusion  collision  diameter  a' 
given  by 


=  1-3  In  1  • 


(62b) 


The  atom-molecule  collision  diameters  wi]l  be  taken  as  the  arithmetic 
average  of  the  atom-atom  and  of  the  molecule-molecule  collision  diameters. 
This  assumption  corresponds  to  the  concept  that  the  collision  diameter  is 
a  measure  of  the  effective  range  of  the  electron  distribution  about  the 
molecule  ox  atom,  and  that  a  collision  occurs  whenever  these  electron 
distribution-  overlap. 

With  the  onset  of  ionization,  several  additional  types  of  collisions 
must  be  considered.  These  are  the  ion-atom,  electron-atom,  and  the  ion- 
ion,  ion-electron,  or  electron-electron  collisions.  The  ion-atom  colli¬ 
sions  may  be  treated  in  the  same  manner  as  the  atom-atom  collisions,  since 
the  spectroscopic  data  for  ionized  molecules  are  available  to  be  used  in 
setting  up  the  Morse  potential  function  between  the  ion  and  the  neutral 
atom.  The  Morse  function  constants  for  these  cases,  which  occur  In  air, 
are  also  listed  in  the  preceding  table.  The  ion-ion,  ion-election,  and 
electron-electron  collision  diameters  all  depend  on  long  range  electro¬ 
static  forces  of  identical  magnitude.  The  criterion  for  the  collision 
diameter  will  be  taken  as  before  ju(cr)  |  =  kT.  The  function  U  in  this 
case  is  just  the  coulomb  potential  between  two  charges,  and  the  collision 
diameter  thus  becomes 
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For  electron-atom  collisions,  the  atom  will  t>e  polarized  by  the 
approach  of  the  electron  and  there  will  result  a  charge-dipole  cype  inter¬ 
action.  In  order  to  calculate  the  magnitude  of  this  interaction,  the 
polarizability  of  the  atom,  a,  will  be  needed.  The  polarizability  is 
defined  as  the  dipole  moment  produced  in  the  atom  by  an  electric  field,  E, 
of  unit  strength.  The  interaction  energy  of  the  induced  dipole  moment 
with  the  field  is  aE2/2,  and  when  the  field  due  to  the  presence  of  the 
electron  is  substituted  for  E,  this  energy  becomes 


U  =  - 


e^a 

2r4 


(64) 


It  will  be  noted  that  the  polarizability  need  not  be  known  with 
great  precision,  since  the  collision  diameter  depends  only  on  the  fourth 
root  of  this  quantity.  Experimental  values  of  polarizability  are  only 
given  for  the  molecular  state  of  oxygen  and  nitrogen.  However,  Joos 
(ref.  22)  gives  an  approximate  method  of  calculating  the  polarizability 
of  atoms  in  alternating  fields  which,  extrapolated  to  steady  state  fields, 
becomes 


where  m  in  this  case  is  the  mass  of  an  electron  and  IQ  is  the 
resonance  potential  for  ionization  of  the  atom.  This  resonance  potential 
is  9.U  e.v.  for  oxygen  atoms  and  10,26  p.v.  for  nitrogen  atoms  (ref.  23). 
Using  these  values,  one  finds  the  polarizability  is  13.2xlC"25cm3  and 
10. 3xl0-25cm3  for  oxygen  and  nitrogen  atoms,  respectively. 

The  collision  diameters  will  enter  the  calculations  to  follow  in 
the  form  of  a  ratio  with  the  collision  diameter  for  two  diatomic  molecules 
at  the  same  temperature.  These  ratios  were  calculated  by  the  methods 
outlined  above  for  those  collisions  which  occur  between  the  major  com¬ 
ponents  in  high-temperature  air.  The  results  are  given  in  table  V  as 
functions  of  temperature.  Wherever  two  or  more  types  of  collision  give 
the  same  result  within  5  percent,  there  is  not  much  point  to  considering 
them  separately,  in  view  of  the  other  approximations  involved.  There¬ 
fore,  for  such  cases,  the  values  have  been  weighted  and  averaged,  and 
only  those  final  averages  are  given  in  table  V.  The  atom-molecule  cross 
sections  S(02-0),  S(N2-N) ,  and  S(H2-0)  are  very  close  and  these  are  all 
designated  by  the  single  symbol  S(N2-N) .  Similarly  the  atom-atom  and 
atom-ion  cross  sections  S(0-0),  S(N-N),  S(N-O),  S(Mf) ,  and  S(0-0+) 
are  all  grouped  under  the  notation  S(ll-H) »  The  atom-electron  cross 
sections  S(N-e)  and  S(O-e)  are  given  the  single  average  value  listed 
under  S(N-e) .  In  the  same  way,  all  the  atom-molecule  and  atom-atom 
diffusion  rrr.  (5  5  euctloiLB  an  given  by  O'  (S2-ft)  a :A  S’ (1-3),  reotp^filvely. 
In  view  of  the  order  of  the  approximations  involved,  the  electron- atom 
cross  sections  for  diffusion  are  not  differentiated  from  the  cross 
sections  for  momentum  and  energy  transport. 
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Viscosity 


The  viscosity  of  air  will  be  calculated  front  the  simple  summation 
formula  for  a  mixture  of  hard  spherical  molecules 

"’i  IPiUiAi  (“) 

i 

In  view  of  the  other  approximations  involved,  a  more  sophisticated 
formulation  than  this  does  not  seem  warranted.  The  constant  5jt/ 32  is 
chosen  to  agree  with  the  value  which,  accounts  for  the  persistence  in 
velocity  of  the  higher  speed  molecules  (ref.  1).  It  will  be  convenient 
to  use,  as  a  reference  value,  the  viscosity  which  air  would  have  at  the 
same  temperature  if  the  molecules  did  not  dissociate  or  ionize 


lo 


5* 

32 


PoHo^o 


(67) 


In  cgs  units,  this  reference  viscosity  becomes 


^o 


1.462x10" 5 


1+112 /T  cm- sec 


(68) 


where  T  is  the  absolute  temperature  in  degrees  Kelvin.  The  ratio  of 
mean  velo^-ties  is 


(69) 


and  the  ratio  of  densities,  where  the  total  number  of  moles  per  unit 
volume  is  the  same,  is 


Pi 


P 


o 


(70) 


Under  these  same  conditions,  that  is,  the  same  mol  concentration,  the 
ratio  of  mean  free  paths  is  (ref.  l) 
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Then  the  ratio  of  viscosity  to  the  reference  viscosity  becomes 


a.  y  Phx.h 

nQ  L  7  Mo  1  \ 


.  Bef°r  Ration  (71)  can  be  applied  to  the  case  where  ionization 
i  Xt  be  obser'/ed  that  tlle  atom-electron  collisions  are  not 

ty  bf  countfd  in  evaluating  the  atom  mean  free  paths,  nor  are  the  ion- 

This  r^becai1!1^  ^  te„Coulltfcd  in  evaluating  the  ion  mean  free  paths.  " 
This  s  because  the  mean  free  path  which  occurs  in  the  derivation  of 

'  i  iS  distajlct  between  those  collisions  which  cause  a 

ei^eJ  i^dl rJSt / Change/n  ffiOIEentum  of  the  particle  being  considered, 

^  1  direction  or  in  magnitude.  Normally  all  collisions  qualify 

in  this  sense,  but  when  a  heavy  particle  like  an  atom  or  an  ion  collides 
with  an  extremely  light  particle  like  an  electron,  there  is  practically 
no  change  in  the  momentum  of  the  heavy  particle  and  this  momentum  is 
carried  intact  to  the  first  collision  with  another  heavy  parti 

A  feature  of  the  ionized  gas  which  should  be  mentioned  is  that  charced 
particles  under  the  influence  of  coulomb  potentials  may  be  deflected  by 
the  cumulative  effect  of  many  long-range  collisions.  Then,  as  pointed^out 
by  Cohen  Spitzer,  and  Routly  (ref.  24),  the  effective  collision  cross 
secexon  is  larger  than  for  a  single  collision  at  short  range.  However 
o^-sMeldliig ^ffect  limits  the  range  of  the  coulomb  potential  and 

ion’  ±TZ»t-i  enSlty+iS  largG  thS  correction  1°  the  collision  cross  sec- 
ion  is  small.  For  the  range  of  variables  treated  here,  collisions  with 

ThusrainP?r  ^eS+Senerally  predomlnate  when  the  charge  density  is  small, 
thus,  in  the  first ^approximation  the  mean  free  path  is  considered  to  be 

ldei.ee  of  single,  independent  collisions  as  in  equation  (71). 

,  .  :Fhf  coefficients  of  viscosity  for  high-temperature  air  have  been 

(68>>  (71>>  Mld  (?2)  the  results  are  i  re- 

!5uatlon  fral1  i  “i  in  "T*  8'  The  terms  ln  the  summation  of 
equation  (78)  will  also  fce  used  in  the  evaluation  of  the  coefficient  of 
thermal  conductivity  which  follows.  nent  01 


Thermal  Conductivity 


A  theory  for  the  thermal  conductivity  of  a  chemically  reacting  gas 
was  perhaps  first  outlined  by  Nemst  (ref.  25).  The  fomn  in  which  this 
theory  has  been  developed  by  Hirschfelder  (ref.  15)  will  be  used  here  to 
ca  culate  the  conductivity  of  high -temperature  air,,  In  this  theory  the 

firsty-na^Tfe^thr°,JSh  the  SaS  iS  treated  ^  two  independent  parts.  The 
first  part  is  the  energy  transferred  by  molecular  collisions,  Ld  this 

mode  of  energy  transfer  is  the  one  responsible  for  the  ordinary  thermal 
conductivity  of  nonreacting  gases.  The  second  part  is  the  en^gy  in¬ 
ferred  by  diffusion  of  the  molecular  species  and  the  reactions  which 
occur  as  the  gas  tends  to  maintain  itself  in  chemical  equilibrium  at  each 


m 
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point.  We  shall  first  turn  attention  to  the  calculation  of  energy 
transferred  by  molecular  collisions.  If  one  invokes  Eucken’s  assumption 
(ref.  26),  that  the  internal  energy  is  distributed  among  the  gas  particles 
independently  of  their  velocity  distribution,  then  the  simple  linearized 
expression  for  the  coefficient  of  thermal  conductivity  in  a  mixed  gas 
becomes 


^  ^  (73) 

where  C-fc  is  the  specific  heat  per  mol  due  to  translational  energy 
^int  4be  specific  heat  due  to  the  internal  energy.  It  is  noted 

that 


ci  -  (Ct  +  Clnt>i 


and 


whence  equation  (73)  becomes 


= 


9 

4  Mi/ 


(7J+) 


(75) 


(76) 


It  will  be  convenient  to  define  a  reference  coefficient  of  thermal 
conductivity  just  as  was  done  for  the  viscosity,  that  is,  the  value  air 

would  have  if  it  did  not  vibrate,  dissociate,  or  ionize.  This  coefficient 
will  be 


or  in  cgs  units 


*o 


19  R_ 

h  Mo  n° 


ko  .  l.36k  %  - 

cm- sec  K 


(77) 


(78) 


where  is  in  gm/cm-sec. 

The  ratio  of  the  thermal  conductivity  coefficient  to  the  reference 
coefficient  becomes 
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The  first  factors  in  each  tern  are  just  the  terms  already  calculated  for 
the  viscosity.  It  may  be  noted  that  the  same  collisions  are  to  be  neg¬ 
lected  in  this  case  as  for  viscosity.  This  is  because  these  electron- 
atom  and  electron-ion  collisions  have  only  a  small  effect  on  the  magnitude 
and  direction  for  the  energy  flux  vector  of  the  heavy  particles.  Actually, 
all  of  the  electron  collisions  could  have  been  neglected  for  the  purpose 
of  computing  viscosity,  for  the  electrons  carry  a  negligible  fraction  of 
the  total  momentum  transferred,  due  to  their  small  mass.  However,  because 
of  its  high  velocity,  an  electron  transports  a  large  chare  of  the  kinetic 
energy,  in  fact  much  more  than  a  heavier  gas  particle.  Therefore  the 
heat  conduction  terms  for  the  electrons  are  needed  in  equation  (79) • 

The  second  mode  of  energy  transfer,  whi ch  takes  place  whenever  the 
gas  undergoes  a  chemical  reaction,  is  due  to  the  diffusion  of  the  chemical 
species.  These  particles  then  react  with  one  another,  giving  off  or 
absorbing  the  heat  of  reaction  and  causing  heat  transfer  which  may  be 
uunaiderably  larger  than  the  ordinary  heat  transier  due  to  molecular 
collisions.  Hirschfelder  (ref.  15)  has  formulated  this  problem,  but  in 
terms  of  the  multicomponent  diffusion  coefficients  which  are  difficult 
to  estimate.  However,  Butler  and  Brokaw  (ref.  2j)  have  shown  how 
Hirschfelder' s  results  can  be  modified  to  make  use  of  the  binary  diffu¬ 
sion  coefficients  instead,  and  their  solutions  are  in  convenient  form 
fur  computatiunal  purposes.  Moreover,  Butler  aid  Brokaw  show  that 
Hirschfelder' s  method  of  predicting  thermal  conductivity  agrees  very  well 
with  the  experimental  results  for  a  number  of  gases  which  dissociate  or 
reeet  chemically  at  ordinary  temperatures.  Therefore  we  shall  be  able, 
tc  apply  this  method  with  some  confidence.  That  part  of  the  coefficient 
of  thermal  conductivity  which  is  due  to  the  chemical  reaction  will  be 
designated  kr,  and  according  to  the  results  of  reference  27,  this  may 
be  expressed  as 


kr 


R 


d  Zn  KpY5 
dT  / 


1 1 31k  (aixJ-aJxi) 

i  j 


(80) 


where  D ' j  is  the  binary  diffusion  coefficient  for  molecules  of  type  i 
into  molecules  of  type  j,  and  the  ap  are  the  stoichiometric  coefficients 
of  components  Ap  in  the  chemical  reaction  written  in  the  form 
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-  0 


(81) 


According  to  reference  1, 


the  binary  diffusion  coefficients  are  given  by 


_i_  _  8  [2  (  No 

nDtJ  3  V  *  )  (RT)i /2  °iJ 


Mat  f  4-  «  — *  —  4*  —  -w  n  „  ,  „ 

w  reference  coenicienr  or  thermal  conductivity  (eq.  (77))  can  be 
put  in  the  form 


95  n/3F  JJ  /pip 
~  N0S0  V  Jfe 


(83) 


one^bta^ns8  eqUati°n  ^  in  ecluation  (80)  and  dividing  by  equation  (83), 


In  order  to  simplify  the  computations  somewhat,  the  differences  In 
mass  between  oxygen  and  nitrogen  atoms  will  be  neglected.  Then  since 
^me,°f^he  collision  cross  sections  may  also  be  given  the  same 'value 
the  double  summations  in  the  denominator  of  equation  (84)  take  on  the 
relatively  simple  form  as  follows: 


For  case  I  (oxygen  dissociation  reaction  only) 


S*(Ng-N)  f [x(0)+2x(Qg)32  x  4x(Wg)|  (  x(N2) 

^3  s0  1  x(0)x(02)  T  x(0)  J  +  ^2  x(02) 


(85) 


for  case  II  (nitrogen  dissociation  reaction  only) 


■  •ifcflftjr  ft-.-  Vgy-t  {-svc.--  •  ^  -I*-. 
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Y  V  »  S>  ^Ng-w)  |  [x(N)+2x(N2)  3^  +  x(0)  |  +  S»  (jf-N)  2x(0) 
i  j  \f3  S0  ^  x(N)x(N2)  x(N2)  J  S0  x(N) 


axid  finally  for  case  III  (ionization  reaction  only) 


IH§- 


i  J 


(N-N)  +  _1_  S»  (N-e~) ~[  [x(N)+x(»+>  ]' 

S°  23°  S  J  x(N)x(N+) 


The  total  coefficient  of  thermal  conductivity  is  just  the  sum,  kn+kr. 
The  values  of  this  sum  were  calculated  for  air  and  the  results  are  pre¬ 
sented  in  table  VI (b)  and  graphed  in  figure  9.  The  coefficient  of  thermal 
conductivity  has  about  the  same  functional  features  with  respect  to  tem¬ 
perature  and  pressure  that  were  observed  for  the  specific  heat  (fig  4) 

The  Prandtl  numbers  were  calculated  from  the  relation 

pr  .  =  1  (52e)  m 

Mk  19  V  R  J  k/ko 

and  they  are  listed  in  table  VI (c)  and  are  graphed  as  functions  of 
temperature  in  figure  10. 

As  noted,  nitrogen  and  oxygen  have  approximately  equal  collision 
cross  sections,  so  that  up  to  the  point  where  ionization  begins,  air  can 
e  treated  essentially  as  a  two-component  mixture  of  atoms  and  molecules. 

L,/®  P°SSil3le  to  characterize  the  system  with  a  single  diffusion 

en"'  ^  the  dlfferential  equations  of  fluid  flow  may  be  analyzed 
■while  keeping  the  terms  describing  heat  transfer  by  chemical  reaction" 
separate  from  the  terms  describing  ordinary  thermal  conduction.  This 
method  is  followed,  for  example,  by  Fay  and  Riddell  (ref.  28)  in  their 
^abf8lSc;±0r,t51e,hfJt  transfer  to  the  stagnation  region  of  blunt  bodies 
“  ti-ight .  In  using  such  an  analysis,  it  must  be  realized 

that  the  thermodynamic  and  transport  coefficients  are  defined  somewhat 
differ entry  than  in  the  present  paper.  The  coefficients  which  appear 
in  reference  28  will  be  called  the  "partial  coefficients"  and  they  will 
be  designated  here  by  a  superscript  prime.  The  partial  specific  heat, 
he  Partial  coefficient  of  thermal  conductivity,  the  partial  Prandtl  num- 
r,  and  the  partied.  Lewis  number,  respectively  are  related  to  the 
quantities  defined  in  this  paper  as  follows  1 


Cp'  =  J^xi(ci+1) 


k*  =  k 


n 


Pr'  =.  n 

Mk' 


(89) 


Le*  = 


DpCr 


Met 


I 

J 


and  ^  e,Ste?  a^dlmensionless^group  **“  *  *"  ^2> 


ZDpR 

koMo 


=  O.309 


S'(N2-N) 


(90) 


The  Lewis  number  is  conveniently  calculated  from  the  expression 


Le,  _  ZDpH  Zty/fc 

2kn/kQ 


(91) 


zr  *  /vTI  thG  calculated  values  for  the  partial  coefficients 

az^t^/trkP^s;in“nrs  s:  md  Le;  « 

not  given  ton  temperature" 

h®r‘  “  '  !ea=\ a  three-component  system  (neutral  atoms.  Ionized  atoms 

ZsTaT^l  £  W^at^Tt1?’  ftl^’  d°  ^  ^  t0 

hand,  decreases  as  dlssoclatLl  proceMHo  S^tlon”^  <V“  ^ 

ss  2:  jet  s-  ^rHTrS,^ 

at  about  the  same  rate  (see  en  (9,9}  <-,■».  fnr\\\  j  j-ncreases 

two  factors  is  ^ota+r  }  q*  (°2J  0r  (90)),  and  the  product  of  these 
iactors  is  relatively  constant:  thus  the  partial  Tcv-ic  , 

approximately  as  M  that  is  i+  *  partial  Lewis  number  varies 

about  1.4  to  ol  '  ’  '  decreases  by  a  factor  of  about  2,  from 

A  list  of  conversion  factors  is  given  in  tabic  vttt 

in  converting  the  tabulated  and  !  1  VI11  for  convenience 

g  tne  tabulated  and  graphed  parameters  into  dimensioned  units 
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DISCUSSION  OF  RESULTS 


It  is  convenient  to  summarize  thp 

the  transport  Properties  of  air  in  «  ^  .  processes  which  influence 

Prandtl  nlbe/vith  A?  wT*1”1  °f  the  varia“ons  in 

pure  diatomic  gas  wlthT  specific  heat  ^  ls  U*e  a 

the  Prandtl  number  for  such  a  mas  is  U/S  or  0  If?,  (88)' 

increases,  vibrational  energy  is  excited  It,  the  1  A,  ”I"erature 

approaches  9/2  R  wh-»le  fv™  th  m°lecules  so  that  cl 

approaches  (23/4)  (En^/M  ^  w  ^  +-  . on  (^6)  the  thermal  conductivity 
18/23  or  0.78.  (  n°/^)‘  F°r  thl"  li£lit  the  prandtl  number  becoLs 

At  still  higher 
temperatures  the  oxygen 
dissociates  and  both 
Up  and  k  go  through 
the  pronounced  maxima 
shown  in  figures  4 
and  9.  To  a  first 
approximation,  k  is 
about  proportional 
to  Cp  just  as  for  a 
nonreacting  gas.  How¬ 
ever,  due  to  the  influ¬ 
ence  of  the  chemical 
reaction,  the  maxim’ im 
in  k  occurs  at 
slightly  lower  temper¬ 
atures  them  the  maximum 
in  specific  heat  as 


.0-f 


Temperoture,  T 

Sketch  (b) 


*  - **'-*-*  w  - \  / 

u  ttt  STM  1,  ,,.1.,, 

tt “'SS  £ ~  L  2  " 

figure  10.  iunction  can  be  seen  in  the  curves  of 

and  nitrogen  dSsocf^Ition^as  n^yet^e^1^^^13  essentiany  complete 
calculated  for  a  two-oo^nt  2/2e  f™’fJh\Pr“dtl  D“ber  caD 
the  mol  fractions  are  x\  ”  s/lST  T?/  con41tl“s  that 

weights  are  related  by  (vh  2  =  ^  ^  that  the  Elolecula-r 

molecules  and  subscript  2  to  ox^gL  atoms^  'f^h1  "Tf*  t0  nitr0B™ 
sections  are  all  assumed  equal/the  j/  L'l/  collision  cross 

Si  une  mean  free  paths  are  related  by 

Xi 


(92) 
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Then  the  Prandtl  number  becomes 


»)©■»-) 


•{2  ^ 


s  0.  72 


which  checks  closely  with  the  values  for  Prandtl  number  shown  in  figure  10 
at  the  transitions  from  oxygen  to  nitrogen  dissociation. 

As  nitrogen  dissociation  proceeds ,  the  Prandtl  number  again  follows 
an  S  shaped  function  of  temperature  for  the  same  reasons  which  were 
outlined  for  oxygen  dissociation  (sketch  (b)).  The  fully  dissociated 
air  is  like  a  pure  monatomic  gas  with  Cp  ‘about  5/2  R  and  k  equal 
15/4  (RTl/fo),  so  that  the  Prandtl  number  approaches  2/3.  The  Prandtl 
numbers  of  figure  10  go  through  this  value  of  2/3  at  a  temperature  where 
the  dissociation  is  essentially  complete  but  the  ionization  is  still 
negligible  (Z  =  2,  fig.  1). 

Up  to  this  point,  the  Prandtl  numbers  lie  within  the  range  from  0.6 
to  1.0  in  agreement  with  the  conclusions  of  reference  29.  However,  when 
ionization  begins,  the  Prandtl  number  may  drop  to  somewhat  lower  values. 
For  very  small  degrees  of  ionization  it  can  be  shown  that  the  ratio  of 
electron  to  atom  mean  free  path  is  about  V ~2  and  that  the  coefficient  of 
viscosity  is  influenced  very  little  by  the  presence  of  the  electrons. 
However,  the  thermal  conductivity  is  greatly  increased  by  the  electrons 
because  of  their  high  thermal  velocity,  and  the  effect  is  illustrated  in 
figure  9  by  the  abrupt  change  in  slope  for  k  where  the  ionization  reac¬ 
tion  first  begins.  If  the  fraction  of  ionization,  e,  is  so  small  that 
chemical  reaction  terms  may  be  neglected 


where  Mx  is  the  atom  mass  and  M2  the  electron  mass.  It  can  be  seen 
that  k  increases  rapidly  with  the  onset  of  ionization  because  the 
factor  (2M1/M2)1/2  is,  large,  about  230.  Accordingly  the  Prandtl  number 
is  found  to  vary  as 


r$\  fr 

KJi  )  \k, 


5  R\  15  R_ 
2  HI/  4  M  :L 
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A  very  small  value  of  e  is  sufficient  to  reduce  the  Prandtl  number 
considerably  below  2/3,  the  value  for  a  pure  monatomic  gas.  As  the 
ionization  proceeds,  however,  the  chemical  reaction  terns  which  were 
omitted  in  equations  (94)  and  (95)  become  predominant  such  that  Cn 
becomes  approximately  proportional  to  k  again,  and  the  Prandtl  nSmber 
levels  out  at  about  0.3  as  shown  in  figure  10.  At  this  point  another 
factor  gradually  predominates,  namely  the  decrease  in  coefficient  of 
viscosity  which  occurs  as  the  mean  free  paths  become  vex-y  short  as  a 
result  of  the  long  range  coulomb  forces  acting  on  the  charged  gas  parti- 
cles.  Consider,  for  example,  the  completely  ionized  gas.  The  mol  frac¬ 
tion  of  ions  (subscript  1)  equals  the  mol  fraction  of  electrons 
(subscript  2)  and  the  collision  cross  sections  are  about  equal.  Then 
the  ratio  of  mean  free  paths  is 


^2 

Ai 


(96) 


where ,  as 


before,  the  electron  collisions  do  not  count  toward  the  ion 

—  IT1L  _  _  1  »  1  /  ■ 


'  - - ■ — kj  ^uuii  1  bUwti 

mean  free  paths.  The  ratio  k/r\  is  approximately  given  by 

k  _  15  _R_  «/2 

t  Mi  Iwi 

and  with  Cp/M  about  5R/M1,  the  Prandtl  number  becomes  very  small 

GP  T]  _  4  1+  V2 

M  k  '  3  V2  V  Mi  =  0,1 


(97) 


.014 


(98) 


as  shown  in  figure  10.  With  the  assumption  that  the  Prandtl  number 
retains  its  usual  significance,  boundary -layer  regions  in  highly  ionized 
flow  should  be  much  better  heat  conductors  than  in  the  fjow  of  neutral 
gas  partides  This  result  is  of  theoretical  interest  even  though  fully 

sS1rS.f60Lfy).a0t  be  “4  at  fllEht  "***«“ 

The  Prandtl  number  functions  of  figure  10  are  similar  to  those 
predicted  by  Kaeppeler  and  Krause  (ref.  30).  Their  calculations  were 
baoed  on  the  low  value  for  nitrogen  dissociation  energy  and  so  the  pres- 
en  estimates  exhibit  more  distinctly  separate  maxima.  This  is  because 
^eparai:ion  between  the  dissociation  reactions  which  occurs  when 
fref^firh^1^  iOV  n±tTO&en  dissociation  energy  is  used.  Greifinger 

air  *+ a  I®*®  S°me  estimates  fo1’  the  transport  properties  of 
air  at  high  temperatures.  He  assumes  constant  collision  cross  sections 

UP  t0  the  P°lat  Vhere  ionization  begins,  his  calcu- 
lated  Prandti  number  is  a  constant  lacking  the  structure  shown  in  fig- 

aie/u*.  however,  reference  31  predicts  the  same  order  of  effect  due  to 
i jnization  which  has  been  discussed  here. 
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Recall  that  the  nitric  oxide  in  air  has  been  neglected  ar.id  that  the 
calculated  themodynamic  properties  of  air  are  not  greatly  affected  by 
this  omission.  It  is  not  obvious  that  tne  transport  properties  of  air 
will  be  similarly  independent  of  NO  formation,  since  NO  is  ionized 
rather  easily  compared  to  N2  and  02  molecules.  Practically  all  of  the 
electrons  which  appear  in  air  at  intermediate  temperatures  come  from 
the  NO  ionization  reaction.  However,  Gilmore  finds  that  the  quantity 
of  such  elections  is  very  small,  the  order  of  0.01  percent  or  less  (ref.  9) , 
and  from  equations  (94)  and  (95)  it  can  be  judged  that  this  amount  is  not 
yet  sufficient  to  alter  greatly  the  transport  properties  of  air. 

At  present  there  a, re  no  direct  measurements  of  the  transport 
properties  for  air  at  the  high  temperatures  considered  here.  One  experi¬ 
mental  comparison  which  can  be  made  at  the  moment  is  to  insert  the  calcu¬ 
lated  coefficients  into  a  theory  for  the  heat  transfer  to  a  simple 
aerodynamic  shape,  and  then  compare  the  results  with  measured  heat- 
transfer  rates.  For  example,  reference  32  compares  the  values  for  heat 
flux  to  the  stagnation  region  of  a  blunt  body  in  high-energy  air  flow 
which  were  measured  by  Rose  and  Riddell  (ref.  33)  with  the  theoretical 
predictions  based  on  the  transport  coefficients  developed  in  the  present 
report.  It  is  found  that  the  calculated  heat  flux  deviates  less  than 
10  percent  from  the  mean  of  the  measured  data.  This  is  a  favorable  result 
in  view  of  the  approximations  involved  in  the  theory  and  of  the  scatter 
in  the  experimental  data  (the  order  of  ±20  percent).  It  should  be  empha¬ 
sized  that  this  comparison  is  not  a  sufficient  test  for  the  reliability 
of  the  calculated  transport  coefficients,  as  it  is  possible  that  errors 
in  the  coefficients  compensate  for  each  other  or  for  approximations  in 
the  heat-transfer  analysis.  Nevertheless,  the  favorable  comparison  is 
reassuring  inasmuch  as  it  means  that  at  least  a  necessary  condition  for 
the  correctness  of  the  theories  is  fulfilled. 


CONCLUDING  REMARKS 


In  conclusion,  the  thermodynamic  properties  (including  compressibility, 
energy,  entropy,  specific,  heat,  and  the  speed  of  sound)  and  the  transport 
properties  (the  coefficients  of  viscosity  and  thermal  conductivity  and 
the  Prandtl  number)  have  been  estimated  for  air  at  high  temperatures. 

These  estimates  were  made  from  approximations  in  closed  form,  and  it  is 
found  that  these  appro ximat ions  give  fairly  accurately  the  thermodynamic 
properties  over  the  range  of  pressures  and  temperatures  for  flight  through 
the  atmosphere.  The  transport  properties  were  calculated  by  methods  which 
have  given  reasonably  good  results  for  gases  at  ordinary  temperatures  and 
for  gases  which  dissociate  and  react  chemically.  Moreover,  the  values 
calculated  for  the  transport  properties  can  be  used  to  predict  the  same 
order  of  heating  rates  that  have  been  measured  at  the  stagnation  region 
of  blunt  bodies  in  high-energy  air  flow.  Therefore  it  is  concluded  that 
the  results  presented  in  this  paper  are  useful  engineering  approximations 
for  the  properties  of  high-temperature  air.  In  particular,  it  is  an 
advantage  that  the  solutions  appear  in  closed  form  because  the  functional 
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relationships  can  be  visualised.  In  addition  the  air  properties  can 
conveniently  be  computed  in  very  small  intervals  to  give  tables  that  can 
be  used  in  obtaining  solutions  1o  real  air  flow  by  the  method  of  char¬ 
acteristics,  or  for  numerical  sclutlcns  of  boundary- layer  flow.  It  is 
expected  that  the  calculated  coefficients  of  viscosity  and  thermal  con¬ 
ductivity  will  prove  to  be  lower  bounds  to  more  precise  calculations  or 
measurements.  The  prediction  that  the  Prandtl  number  for  highly  ionized 
flow  will  be  small  compared  to  unity  is  of  interest,  for  it  implies  that 
boundary- layer  regions  in  such  a  gas  will  be  very  transparent  to  heat 
flux. 


Ames  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Moffett  Field,  Calif.,  Nov.  18,  1957 
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TABLE  I.-  PARTITION  FUNCTION  CONSTANTS  FOR  THE  MAJOR  COMPONENTS  OF  AIR 
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S(N2-N) Is(N-N) (S(N-e) Is(e-e) |S* (N2-N) 


So  I  So 


500 

1,000 

1.500 
2,000 

2.500 
3,000 

3.500 
4,000 

4.500 
5,000 

5.500 
6,000 

6.500 
7,000 

7.500 
8,000 

8.500 
9,000 

9.500 
10,000 

10.500 
11,000 

11.500 
12,000 

12.500 
13',000 
13,500 

14,000 
14' 500 
15,000 


38.4 

34.9 

33. 7 

33.2 

32.8 

32.6 

32.4 

32.3 
32.2 
32.1 
32.0 
32.0 

31.9 

31.9 

31.9 

31.8 

31.8 

31.8 

31.8 

31.8 

31.7 
31-7 
31-7 

31.7 

31.7 

31.7 
31-7 

31.7 
31.6 
31.6 


0.946 

.920 

.889 

.886 

.846 

.830 

.815 

.803 

•  792 

•  782 

•  773 
.764 

•  757 

•  750 
.743 
.737 

•  731 
-725 
.720 
.715 
.710 
.706 
.701 
-697 
.693 
.689 


0.894 

.838 

.785 

.742 

.705 

.676 

.650 

.628 

.608 

.591 

.575 

.561 

•  548 

•  536 
.524 
.514 
.504 
.495 
.486 
.478 
.470 
.463 
.456 
.448 
.443 
.437 
.431 
.426 
.420 
„4l5 


0.397 

•  380 
.366 

•  353 
.342 

•  331 
.321 
.313 
.304 
.297 
.290 
.283 

I  .281 

!  .270 
.266' 
.261 
.256 
.252 
.247 
.243 


89.9 

75.6 

64.5 

55.7 

48.6 

42.8 

37.9 

33.8 

30.4 

27.4 
24.9" 

22.7 

20.8 
19.09 
17.60 
16.27 

15.10 
14.04 

13.09 

12.24 


0.877 
.843 
.817 
.794 
.775 
.759 
.745 
.733 
,722 
.712 
.703 
.695 
.688 
.681 
.674 
.668 
.662 
.657 
.652 
.647 
.642 
.637 
•  633 

.629 

.625 

.621 


0.761 

.703 

.652 

.611 

.578 

.551 

.527 

.507 

.489 

.473 

.458 

,445 

.433 

.422 

.412 

.402 

.393 

.385 

.377 

.370 

.363 

.356 

•  350 
.342 

•  338 
.332 

•  327 
.322 
.316 
.312 


S(Na-N)  =  S(02-0)  =  S(Na-0) 

S(K-N)  =  S(0-0)  =  S(N-O)  =  SCN-FM  =  S(0-0+) 
S(N-e-)  =  S(O-e-) 

S*(N2-N)  -S*(02-0)  =S»(N2-0) 

S*(N-N)  =  S*(0-0)  =  S'(U-O)  =  S’  (K-ir)  =  S»  (0-0+) 
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(a)  Coefficient  of  viscosity 


Ratio  rj/T)c 


Reference 
coefficient, 
tj0  ,  in  10 "6 


100  10 


Pressure,  atmos’’ieres 

|  | lb  sec  gm 

1.0  0.1  0.01  0.001  0.0001  ft2  cm-soc 


500 

1,000 


1.500 

2,000 

2 . 500 

3,000 

3.500 

4,ooo 

4,50*0 

5,000 j 1.022 


1.000  1 1.000 


1.000 


1.001 


10,000 
lv ,  500 )  1 ,  4lo  j  1 . 457 
1] ,OOo|l.467|1.515 
11,500 
12,000 
12,500 
13,000 
13,50011.61711.528 
14 , 000 1 1 . 628 1 1 . 508 
14,500 
15,000 


1.006 

1.020 

1.033 

1.051 

1.086 

1.148 

1.229 

1.294 

1.332 

1.371 

1.392 

1.405 

1.404 

1.309 


1.352 


1.542  1. 


1.287 
1.190 
I.065 
•  908 
.741 
.580 
,443 
.306 
.254 


1.010 

1.022 

1.038 

1.074 

1.146 

1.228 

1.276 

1.317 

1.342 

1.355 

1.354 

1.328 

1.267 

1.162 

1.004 

,806 

.597 

.415 
•  279 
.1918 
.1363 
.1061 
.0843 
.0615 


1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.010 

1.024 

1.055 

1.128 

1.209 

1.257 

1.289 

1.309 

1.315 

1.291 

1.220 

1.085 

.871 

.615 

.384 

.226 

.1362 

.0930 

.0688 

.0598 

.0569 

.0572 

.0576 

.0598 


1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.011 

1.032 

1.096 

1.181 

1.227 

1.259 

1.277 

1.273 


0.558 

.868 

1.100 

1.293 

1.461 

1.612 

1.751 

1.879 

1.999 

2.11 

2.22 

2.32 

2.42 


i=-' 

fc_  m  Jt— 


;  1.086 
.841 
.534 
.280 


.1397 
.6788 
.0556 
.0466 
.0435 
.0436 
•  0453 
.0471 
.0506 
.0543 
.0581 


2.61 

2.69 

2.78 

2.86 

2.94 

3.02 

3.16 

3.17 

3-24 

3.31 

3-38 

3.45 


■3  59 


3-58 

3.65 

3.71 


1011 
1062 
1112 
1159 
12u4 
1247 
1289 
1330 
1370 
1408 
1446 
1462 
1518 
1552 
1586 
1620 
1652 
168  4 
1716 
1747 
1777 
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(b)  Coefficient  of  thermal  conductivity 


Ratio  k/k0 


Pressure,  atmospheres 


0.01  0.001  jo.  0001 


500 

1,000 

1,500 

2,000 

2,500 

3,000 

3.500 

j),000 

4.500 

5,000 

5.500 
6,000 

6.500 


12.000 


1.021 

1.100 

1.150 

1.177 

1.256 

1.421 

1.941 

2.69 
3.22 
3.07 

2.46 

1.930 

3.35 

4.69 
6.31 
8.21 

9.86 
10.90 
10.88 

9.87 

8.33 

6.84 

5.59 

4.79 

3.34 
11.11 
17.62 
20.0 

23.3 

25.3 


1.021 

1.100 

1.150 

1.177 

1.317 

1.928 

3.15 

3.94 

3.06 

1.997 

2.91 

4.53 

6.98 

9.97 


12.48  ]12.24 
13.19  |  7. 80 


11.55 

8.79 

6.38 


5.10 

3.26 

10.09 


4.07  15.44 

Hi  15  19-14 


1.021 

1.100 

1.150 

1.251 

2.50 

5.48 

3.96 

1.600 

3.32 
7. 18 

13.71 

18.74 

15.39 

8.32 
5.92 

3.42 

11.99 

16.83 

23.6 

32.0 

41  iO 


1.021 

1, 


1. 

1. 

4. 

5. 

1. 

2. 

7. 

16.63  125.8 


13.44  25.1 
16.25  31.2 


19.56  38.2 
23.0  45.0 


27.0 

31.3 
35.7 

42.1 

45.3 


52.2 


62.6 

65.3 

64.2 


52.9 

63.2 

70.9 

73.2 
69.8 
61.5 
51.1 

38.7 


33.1 


1 

1 

1 

3 

5 

1 

3 

15. 

30. 

11. 


3.54 
7.30 

14.69 

29.0 

54.3 

87-6 

110.9 

97-5" 

59.1 

30.2 
16.54 
11.41 

9.54 


8.98 

9.08 

9.50 

9-98 

10.73 

11.51 

12.31 


Reference 
coefficient, 
ko,  in  10“6 


Btu  |  watt 


ft  sec  °Rl cm  °K 


5.84 

9.10 

11.53 


*T  o  c;c; 


15.31 

16.90 

18.35 

19.69 

21.0 

22.1 

23.3 

24.3 

25.4 

26.4 

27.3 

28.2 

29.1 
30.0 
30.8 

31.6 

32.4 

33.2 
34.0 

34.7 
35-5 

36.2 


954 

1053 

1143 

1227 

1305 

1379 

1449 

1516 

15.30 

1642 

1701 

1759 

1814 

1868 

1921 

1972 

2C2C 


36.9 
37-6 
38.2 

38.9 


2070 

2120 

2160 

2210 

2250 

2300 

2340 

2380 

2420 


'••■Sr  ' 


mm 


52 
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TABLE  VI.-  TRANSPORT  PROPERTIES  OF  AIR  -  Concluded 


(c)  Prandtl  number  1 

_ _ I 


T, 

Pressure,  atmospheres 

°K 

100 

10 

1.0 

0.1 

0.01 

0.001 

0.0001 

500 

0.738 

0.738 

0.738 

0.738 

0.738 

0.738 

0.738 

1,000 

.756 

•  756 

•  756 

.756 

•  756 

.756 

.756 

1,500 

.767 

.767 

.767 

.767 

.767 

.767 

.767 

2,000 

.773 

,773 

.773 

.166 

.724 

.668 

.614 

2,500 

.762 

•  751 

.696 

.645 

.611 

.  654 

'  .771 

3,000 

.740 

.680 

.627 

.636 

.740 

.745 

•  714 

3,500 

.678 

.631 

.6  60 

.744 

•  737 

•  658 

.606 

4,000 

.640 

.662 

.762 

•  759 

.619 

•  580 

•  587 

4,500 

.654 

.743 

.752 

.610 

.578 

.611 

.764 

5,000 

.702 

.767 

..611 

.581 

.624 

. 799 

•  993 

5,500 

.748 

.620 

•  583 

.617 

.785 

.989 

.871 

6,000 

.763 

•  592 

.602 

.736 

.969 

.891 

•  384 

6,500 

.610 

.592 

.673 

.906 

•  955  "] 

•  363 

.348 

7,000 

•  593 

.620 

.796 

.986 

.830 

.346 

.337 

7,500 

•  595 

.688 

.927 

.969 

.350 

.334 

•  330 

8,000 

.620 

.788 

.983 

.648 

.332 

.328 

.316 

8,500 

.666 

.891 

.943 

•  335 

.324 

.321 

.276 

7,000 

•  730 

.961 

.SO7 

^321 

*307 

.1987 

9,500 

.806 

.  966 

•  330 

.314 

.31 6 

.273 

.1140 

10,000 

.886 

.872 

.308 

.310 

•  313 

.210 

.0577 

10,500 

•  937 

.310 

.301 

.309 

.284 

.1427 

.0312 

11,000 

.955 

.294 

.296 

.303 

.246 

.0870 

.0207 

11,500 

•  947 

.284 

.295 

.293 

.1945 

.0503 

.0157 

12,000 

•  908 

.277 

.293 

,276 

.1409 

.0321 

.0132 

12,500 

.728 

.272 

.290 

.250 

.0949 

.0213 

.0120 

13,000 

.275 

.272 

.284 

.215 

.0634 

.0166 

.0115 

13,500 

•  251 

.270 

.276 

•  1733 

.0416 

.0142 

.0109 

14,000 

.245 

.269 

.263 

.1338 

.0293 

.0130 

.0109 

14,500 

.241 

.265 

.237 

.0903 

.0202 

.0119 

.0109 

15,000 

.238 

.263 

.220 

.0719 

.0119 

.0114 

.0109 

0.01  0.001  0.0001  100 
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Domains  of  pressure  and  temperature  for  the  major  chemical  reactions 


Figure  7*-  Stagnation  temperature  and  pressure  in  air  as  a  function  of 

altitude  and  velocity. 


function  of  temperature. 
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